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Supplementary notes 

This report consists of the results of research work carried out «t 
the Zentralstelle fUr Geo-Photogrammetrie und Fernsrkundung der 
Deutschan Forschungsgemeinschaf t (Central Laboratory of geophoto- 
gramme try and remote sensing of the German Research Council) 

Summary:, 

Based on LANDSAT-l data (photographic products and CCTs) methodological 
and application oriented investigations vere carried out. The geological 
Interpretation on data exhibiting the Italian peninsula ^lineament eva- 
luation) led to the recognition of tectonic features which are explained 
by a clockwise rotation of various blocks along left-handed transform 
faults. These faults can be interpreted as resulting from shear due to 
main stress directed north-eastwards. 

A landuse map of the raountainuous regions of Italy was produced on a 
scale of l;25o.ooo. Several levels, I ,II and III categories could be 
mapped . 

For the digital treatment of MSS-CCTs an image processing software was 
written in FORTRAN tv. The software package includes descriptive 

statistics and also classification algorithms. Statistical investigations 
on redundancy and user-oriented apolication of LANDSAT-data led to inves- 
tigations concerning data reduction and to analogue-digital data pro- 
cessing techniques. At least laser— optical processing techniques fFourier- 
ana lysis) applied on LANDSAT-data were considered. 
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PART I 

SCIENTIFIC AND APPLICATION-ORIENTED INVESTIGATIONS 
■| (GEOLOGY, TECTONICS, LANDUSE-MAPPING) 
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INTRODUCTION 

This final report on LANDSAT-1 investigations is the re- 
sult 03? activities conducted during the years 1973 - 1975. 
With respect to our research activities in the LANDSAT 
program and in the field of remote sensing in general, 
i this report must he considered an intermediate work state- 
ment of progress to date. 

It is necessary to continue and to intensify future v;ork 
on the following aspects: methological investigations 
concerning the applicability of LANDSAT data for various 
earth scientific disciplines; development of data handling 
techniques, such as digital evaluation of LANDSAT data; 
and the analog electronic treatment of images. From a p-ure 
scientific point of vbw, LANDSAT data have delivered a 
great deal of new information. 

With regard to LANDSAT-2, we are optimistic that further 
scientific results can be derived by evaluation data 

I 

acquired over a longer period of time. 

Pure scientific investigations -lead automatically to 
practical applications and the practical use especially 
of LANDSAT data. In order to define and to identify poten- 
tial applications of spaceborne remote sensing technio_ues, 
a first-order task is to bring to the attention of poten- 
tial users the importance of such data. 

The Zentralstelle ftir Geo-Photogrammetrie und Fernerkundung 
has been foimded and is funded by the German Research Council 
in order to promote remote sensing activities in Germany and 
to help establish this new technology for various potential 
applications. The enclosed book on ERTS and partly on SKYLAB 
(" Die dritte Ent eckung der Erde ”) must be considered as 
a step tovrard achieving this goal. The distribution of this 
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book has already resulted in increased interest by potential 
\isera, such as ministries and other agencies involved in 
emrironmental planning and surveillance. 

The LANDSAT-1 data provides us, for the first time, with 
very high resolution sequential multispectral satellite 
data. Our first approach in the evaluation of this data 
has been carried out xmder earth scientific aspects using 
more or less conventional interpretation techniques. These - 
activities consist of structure-tectonical investigations 
and also the evaluation of LAMDSAT data for land-use 
' planning purposes. 

Investigations into the more kethodological aspects con- 
cerning the digital treatment of LAUDSAT data has been 
intensified after receiving the first tapes near the end 
of 1973 . . 


Within the scope of the above-mentioned- investigations , 
the following activities were conducted: software develop- 
ment for the treatment of MSS-data, the application of 
analog digital and optical techniques. 

According to our various activities, the follov/ing reporo 
is subdivided into two main parts entitled: Part I, 
scientific and application-oriented investigations and 
Part II, methodological investigations. 
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A) GEOLOC-IC EVALUATION OF CENTEIAL ITALY 


The geological investigations described are based on 
conventional photo interpretation techniques applied to 
LANjDSAT imagery. For some demonstration pxirposes also 
jE^EP data v/ere applied. 

The, methods described were used to investigate the tectonic 
s'ignificance of identifiable linear features by utilising 
multistage imagery and the results of fieldvrork (Fig- 1). 
The geological field investigations are restricted to some 
areas near Volterra and Siena in the Central Tuscany. 


AREA 



ALPS. 

APENNINES 


TUSCAN-UMESIAN 

APENNINE 


TUSCANY 
N of Vol terra 
of Siena 


Fig, 1 Scheme for the geologic-tectonic research in the 
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I '1 • Methodological Aspects of Space-Imagery Interpretation 

,- 1^1 Multitime Image Interpretation 

I : 

•The ability to recognise and interpret linear features 
within mountain regions strongly depends on shadow effects 
provided by the distinct relief. As shadow effects are 
also a function of Sun angle. ERTS images offer optimum 
possibilities for the statistical comparison of lineament 

maps with respect to seasonal illumination conditions. 

• 

Figure 2 shovrs a statistical representation of lineament 
directions derived from three sequential ERTS images of 
the Molise, northwest of the Gargano Peninsula, Italy 
(ERTS-1 , 8.8.72, 27.8.72, 5.2.73). Comparison of the three 
histograms shows that all possible linear-feature directions 
a-rej present in each ERTS frame. The distributions of the 
maxima, however, have shifted to a certain degree. Depending 
on the azimuth and elevation of the Sim, this effect has to 
be taken into consideration, especially when performing . 

I tectonic interpretations on the basis of statistical linea- 
ment distributions. 

■In order to minimise the above effects, it is necessary to 
superimpose the results of as many multitime lineament 
• evaluations as possible as in the following table. 

' I ■ - - 

•The: significance of the distribution maxima obtained by 
combining sequential data is discussed in Chapter 4. 
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"Fig. 2 Histograms of lineament distributions derived from 
sequential images 

! 

The need for multitime evaluations is also apparent from 
Figure 3. Many linear features could be identified and 
traced over large areas by using images from different 
seasons. The reason for this is that, in Central Europe 
at least the detection of lineaments depends on morpho- 
logy and vegetation. When using multitime imagery, 
regional patterns become more obvious. 
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Flg«3 Composed lineament map of ■temporal ■ imagery 
LANDSAT-1 imases 
MOUSE/ ITALY 

Li neat Ions (MSS 7) 

9, Aug. 1972 . 

27. Aug. 1972 

5, Feb. 1973 

— —overlapping llneations 
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1,2 

By using morphological featvires for the evaluation of 
synoptic satellite imagery , it is possible to obtain 
information on tectonic movements along fracture systems. 
Regions comparable in terms of petrography and tectonic 
•s'tress fields are characterised by a similar morphological 
t|ai|.tern if the same climatologbal conditions can be ass\maed. 
If a displacement of such comparative units is detectable 
along a linear feature in a first approach, it can be des- 
, cribed as a fault. A more detailed description of the fault 

i * , 

is only feasible if the dipping direction of the strata 
can be detected, because vertical displacement of tilted 
’ strata results in a pseudo-horizontal displacement at the 
surface. Vertical movements can be determined only for 

. .fold stru ctures, basin or graben fringes. Sigmoi dal 
. oriented valleys and mountain ridges can mask younger 

^horizontal movements. Furthermore, it can be observed 
that fold structures cut by faults often show displace- 
p'ents with respect to the movement of- the fault. 

Examples of the above phenomena are sho^vn in Figure 4^, 
which represents part of a Skylab image (S 19 oA, Septem- 
j' ber 1973) shov^ing the Umbrian Apennines betvjeen the 

‘ Lago di Trasimeno and the Adriatic coast. 

! i 

.j Figure 5 clearly shows the feasibility of thematic mapping 
of geologic tectonic phenomena. This map, which covers the 
the Latium volcanic region and the northeastern area bet- 
ween the Umbrian Apennines and the Abruzzi, shows the 
advantages that synoptic satellite imagery offers in this 
field. Besides different volcanic phenomena and the strike 
and dip of strata, regional relationships between distict, 
identifiable linear features are also evident. 

Two types of faults in particular can be differentiated: 
transform or shear faults and vertical fractures. 

In this area, the shear faults from a narrow-angle system. 
The first, approximately north-south oriented,, is represen- 
ted by a zone which cuts the Apennines east of the well- 
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knovn Rieti line. Along this fault, horizontal movenent 
can be detected, viith a relative dislocation of the western 
part towards the south. Forming an angle of about 45° t this 
system is cut by southwest-northeast directed shear faults. 
For the north-west oriented faul system, clockwise rota- 


tion of the western and eastern blocks is assumed. 



Fig. 4 Umbrian Arc, Skylab S 190 A; Examples for de- 
tecting faults by morphological features 
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Fvirther detectatle nortlivest- southeast striking vertical 
fault systems are correlated with the chain, of Tuscan- 
Latium volcanoes. These vertical faults, together with 
the above northeast-southwest oriented transform faults, 
form a system of be and ac joints that were originated 
by a southwest-northeast directed stress field. 



5 ; 


r:) 





Geologic tectonic map of LATIUM, ITALY 
LANn$AT-l, 6. Feb. 1973, MSS 7 
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• Lineament Pat'terns of the North and Central Apennines 


Owing to the synoptic nature of satellite data, the exten- 
sive major fault systems could he identified ‘and explained 
in relation to stress directions and fold systems. 



'Fig. 6 Main faults in Italy revealed by ERTS-1 
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Regional Lingamgnts gyal'^ted from 

I 

Figure 5 shows a view of the main lineations in the 
northern and central Apennines* The arrows indicate the 
senses of the horizontal displacements (see Chapter 1) . 

The statistical distribution of the transform and vertical 
faults corresponds with the- results derived bv aerial 
photo-interpretation and fieldwork (Fig. 7 ). Of importance 
for the further interpretation are the distribution peaks 
at 4o^, 7o°, 13o° and l6o° (see Chapter 4). 

Intensive ground studies have proved that these directions 
can be correlated directly with an older tectonic system, 
\rith a~axis 7o° and b-axis l6o°, which was rotated clock- 
wise from a primary system (b l4o°), and with a younger 
system again in the primary stress with a = 4o*^ and b -■ 
l4o°. 

The excellent correlations betvreen the results derived 
from satellite images, aerial photographs and ground tr-uth 
serve to demonstrate the feasibility of the evaluation 
methods applied. 

MIHL PSDT0GRRPH5 i ERTS-l-DflTB 


Fig. 7 Rose-diagram of fault 
directions revealed by 
aerial photographs and 
ERTS-1 data 
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Apennines 

The complex tectonics of the Umbrian arc are demonstrated 
by Fig. S 'which shovrs the main lineaments and accompanying 
minor lineations. Taking into account also the information 
derived from ERTS data (Fig. 7), the same tectonic systems 
a.s above (a.^ = 7o°, b.j = iSo^^j a£ = 132 “ 1^°^) could 

be derived, In addition to these main lineations, shear 
faults at 1o® and llo*^ can be observed. 


E 12 -/oD 



Fig. 8 Photo-tectonic interpretation, Tuscan -Umbrian 
Apenines 
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Some i.nformation could also be derived concerning the age 

*of the various systems. A sequential description could be 
dro-wn up for the latest movements. A model for a ''tectonic 
schedule” is illustrated in Figure 1o and can be described 
as f ollov/s : 

1) Beginning of the folding phase, which is characterised 

northv/est- southeast oriented fold axis and a per- 
pendicular striking ac-system (primary l4o° b and b^). 

2) ¥ith respect to the fiurther compression, displacement 
of the various blocks occurs along the ac-s 3 '‘stem. 

Hiis results in a relative displacement tovrards the 
northeast of the northwestern part. 



Fig, 9 Fracture development in the Tuscan Umbrian Apenine 
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3) Besides north-south oriented transform faults, a con- 
tinuation of the shearing process causes a north east 
directed displacement of the southeastern part. The 
repeated displacements of the folds and overthrusts 
from the "Umbrian Arc". 

4) In a last phase, vertical movements occur especially 
along -70*^ , 11 o° and l6o° striking planes. Along these 
fractures, the region of the Latium and Tuscan vol- 
canoes was depressed. 

General Tectonic Model for the Northern and Central 
3.1 Apennnlnes 


As is clearly shown in Figure 7, the directions of the 
main lineations of the region between the Arno R.iver and 
the Abruzzi Mountains differ from those in the northern 
and southern areas. The tectonic unit of the Tuscan-Umbrian 
Apennines is- bordered by the follov;ing fracture zones; 

- In the north, a transform fault striking about 8o° and 
following the Arno Valley. Along this fault, the 
Tuscan-Umbrian block was displaced to the east. 

- In the east, a north- south striking transform fault, 
beginning east of the city of Latina and ending at 
the Adriatic coast northwest of Ancona. The relative 
movement of the Tuscan-Umbrian -unit was directed to 
the south. 

A cloclntf-lse rotation of the Tuscan-Umbrian Apennines 
results from the combination of movements described. 
Diagrams obtained from fieldwork show tv/o rectangular 
joint systems: and older (b.j = l6o°) and a younger 
(b2 = l4o°) system. In his theory of the orogenesis of 
• the Apennnines, Wmderlich (1966) postulates an anti- 
clockwise rotation of the main stress from west to 
southwest. Such a rotation could have caused the diffe- 
rence between the older and youiger b-axes. As regards 

-JS- 

see next chapter 
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tectonic moveraents revealed by satellite data, it must, 
ho-wever, be pointed out that such studies mainly high- 
light features of relative movements. 


Assuming the Tuscan-Umbrian region to be a static block, 
the forming stress v/ould have shifted from west to south- 
Iwest. However, taking the direction of stress as a stable 
eli^ment, the Tuscan-Umbrian block must be rotated clock- 
wise (Fig. li). This second possibility is confirmed by: 


- the movements along the bordering main transform faults, 

- the internal fault pattern of the block, 

I 

- the overthrusts in the Abruzzi Mountains, 

( I ’ 


I , ’ 

- the reduction in overthrusting intensity from north to 



Fig. lo Tectonic model of xthe.;Tuscan - Umbrian Apennines 



BLATT 16 


ZENTRALSTELltE FOR GEO -PHOTOGRAM METRIE UNO -FERN ERKUNDUNG- 
^ I 

' i ! I 

r ’ 


3 Centra]L_Apennines_ (Abruzzi^M^ 

r\ — — « 

^ ^he tectonic concept, described for the Tuscan-Umbrian 
j Apennines based on the relative movement observable in 
! ERTS and Skylab data has to be seen in connection with 
the neighboinring areas north and south. 

Based on the tectonic map (Fig. 7)* a movement of the 
Central Apennines, the Abruzzi Mountains, similar to 
that of the Tuscan-Umbrian block can easily be derived. 
But the Abruzzi are further southwest than the Tuscan- 
Umbrian Apennines. This could be caused by a combination 
of clockwise rotation and northeastward drift tendency of 
the Tuscan-Umbrian Apennines, whereas the Abruzzi block 
was a relatively stable element. 


3 . 3 

I 

In the west, the Northern Apennines are cut off by a NNE 
striking transform fault beginning in the Gulf of Genoa 
and forming the eastern border of the Monferrato. The 
southern limit is the Arno Valley fault, which borders 
the Tuscan-Umbrian Apennines. 

According to the relative movements along the two bordering 
transform faults, the Northern Apennines must also have 
been rotated clockwise. In this case too, the rotation 

dould be combined with drift in a northeastern direction. 

: I 

These movements may explain the triangular shape of the Po 
' ! 

Valle3’’, the northwestern part of the Northern Apennines 
tleing moved further northeast than their southwestern part. 

A sketch of above tectonic model, evaluated from ERTS and 
Skylab data, is shown in Figure 12. The "block-bordering” 
transform faults can be interpreted as a shearing system 
with a main stress lying S¥-NE. This stress caused a general 
drifting of the Apennines to the northeast. The drift 
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intensity decreased from north-west to southeast and caused 
clockwise rotation of the three blocks of the Northern 
Apennines, the Tuscan-Umbrian Apennines, and the Abruzzi 
Mountains. 

I 

I 



Fig. 11 Main fault systems in Central Italy 
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‘4^ Re latiionsliiT) betv;een Northern Apennines Po Valiev and 

i " * " ' • ' 1 

' * |the Southern Alps 


It is understood that the extensive tectonic fractures must 
have their origin in the deeper crystalline basement. These 
large scale movements result in a system of fractures re- 
vealed only by a subtitle lineament pattern in the sedimen- 
tary cover, interpretable on the basis of space and aerial 
images. 

.However, the exact location of the fracture planes, revealed 
by the lineaments, is difficult to detect. 


, But for the geologist, the knowledge of the relative position 
of the fracture planes, obtainable through evaluation of 
jlbiotollneaments is, especially, of great importance. 

I |- 

Iii the area of the northern Apennines, SW-NE and M-SE 
striking parallel fractures are known and can be traced as 
. lineations through the bordering Po Valley. Some of these 
j lineations even extend into the southern Alpine body. 

> Fhom geological investigations, it is known that swales 
‘ and ridges existed during Triassic time. These features 
resulted from KNE-SSW striking fractures in the crystalline 
■ b-Asement (Auboin, 1967; De Jong, 1967; Nithack, 1968). 

Furthermore, it is proven that there is a stepwise downward 
j progression from the southern Alps to the Po Valley, created 
’ by V7NV/-ESE stepwise faults. 


J 


i'! 


! I 


i' I 


.•t can be expected that the above-mentioned prevalent verti- 
cajl, movements within the crystalline basement are present 
also in the basement of the Po Valley area. 

This conclusion is reached on the basis of the lineament 
evaluation (Figure 15’) and the morphology of the Quaternary 

map of Figure 12 (Beder* ke 


basin, described by the 


Wunderlich, 1968). 

The characteristic morphology of the Quaternary’' basin in 
the region of the ^o Valley suggests that the large scale 
lineaments exist throughout the Apennine-Po Valley and 
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J southern Alps and that they are an expression of actual 
' |£aults. 

Based on this conclusion, a schematic block diagram could 
be constructed for the above-mentioned area(Fig.1'4 ). 
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f 5, Coimnarison of Lineament Pattern with Tectonic. 

j ■' 

Ground Truth Measurement 


Satellite imagery can be regarded as an optimum base for 
a large scale lineament and structure analysis.. By 
utilizing the unique synoptic viev/ over large areas, the 
interpreter is able to delineate linear features and to 
recognize subtle structural elements on the earth sur- 
face, detectable by aircrafts or by ground investigations. 

* \ 

Because of the above-mentioned advantages of space imagery, 
it is not surprising that many scientists, and specially 
geologists, have concentrated their efforts on preparing 
detailed lineament maps. This kind of evaluation is of 
great scientific importance, especially vrith respect to 
the understanding of inter-regional structural relation- 
ships. In addition to the above-mentioned aspects, a 
further evaluation procedvire is to verify the interpreted 
structural elements from space imagery by intensive neld 
work. Within this approach a first order task is to assess 
the influence of interregional structures on .the meso- 
tectonic fabric. 

* ► 

The mesotectoni'c fabric of rocks is the result of certain 
tectonic stress conditions, D^iring geological time, stress 
patterns have changed their orientation and intensities 
several times. The above-mentioned changed stress conditions 
often are reflected by typical symmetrical configurations 
of joints. A careful interpretation of the fabric enables 
the geologist to draw conclusions on former stress condi- 
tions. That is not only of scientific, but also of great 
economic importance. Quite often a close relationship 
exists between a specific tectonic situation and the occurrence 
of ore deposits. In order to investigate the possible influence 
of lineament features to be seen on LANDSAT dat^to the 
raeso-fabric, a test site in the southern part of the 
TUSCAITY near SIENA was chosen. 




rvp THB 
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Fig. 15 t 

Geological Map of the test 
surroundings 

Signatures in darker grey; 



rea and its 

Mesozoic rocks 
Oligocene rocks 
Mlocene/pllocene rocks 
Faults 
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5 . 1 Geolog2_aJ .ihe_Test_Slte 

The groimd truth area is part of a very important 
orogenetic structure with Apenninic strike direction, 
whose nature is still open to discussion. This struc- 
ture is known as Cetona - Orsaro - line ( the most 
suitable name without genetic significance ). 
Morphologically, the Cetona-Orsaro-line forms a 26o kra 
long ridge of variable altitude and extension. In the 
North of the Arno River, it represents the main crest 
of the Apennines. With respect to geology, the ridge is 

I 

built up in its southern part by Mesozoic limestones, 
dolostones and alternating cherty— calcareous-marly strata 
and in Its central and northern parts mainly by the 
Macigno sandstones of Oligocene age. ( Fig. 15 ) In 
general, the structure in the mapping area could .be ex- 
plained as a sequence of two overturned anticlines being 
_ replaced by thrust faults. The original fold-structures 
appear now to be disturbed by synlcinematic shearing and 
postkinematic block -faulting. Such block-faulting revived 
the old fracture pattern and created the horst-like up- 
lifting of the ridge as a whole. The ridge is divided by 
several horizontal faults transverse to the strike of the 
folds C Fig. 16 ). The largest displacements have .taken 
place on ac-direction, as revealed by satellite images 
and aerial photographs ( ’Fig.16.,17) . On the contrary, 

, shear directions are less important with respect to large 
scale displacements. 
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Fig. 18 

I<ANDSAT image of Tuscany 

- : test area rj.;>i:OI>IJOG;CLTTV OP THE! 

PAGE! TO POOH 
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5 j. 2 Evaluation of a LAM)SAT-1 Image of the Tuscany^ with 
Special Regard to the Gro\md Truth Area 

The following Figure 17 is based on a LAM)SAT-1 image of 
the Tuscany taken on August 12, 1972 ( Fig. 18 ). The 
ground truth area is demarcated by dotted lines. The 
interpretation map ( Fig, 17 ) shows all the large 
lineaments touching the calcareous ridge situated between 
Rapolano and Montepulciano. The important ESE-‘WM^’7- striking 
horizontal faults can be followed almost to the Th 3 nr’renian 
coast. In the same way, the bc-fault~ traces ( representing 
reverse faults) reach far to the North and South. In aerial 
photographs ( Fig. 16), a third important fault direction 
striking II 0 ® can be recognized which does not appear on 
the LAM)SAT image. This is probably explained by the fact 
that scan-lines run nearly parallel to this direction. For 
, this reason the recognition of those linears is impossible. 

In addition, the LANDSAT image indicates some lineaments, 
unrecognizable by ground truth or by evaluation of aerial 
photographs. Examples are the lineament ¥ of the ridge, the 
result of the uplifting of the ridge, and the other one N of 
Tre^wanda, the resvilt of the displacement of the ridge by 
about 4oo m transverse to its strike direction. Another one, 
only in part recognizable on aerial photographs, but well- 
identified in the satellite image, marks the boundary of the 
V/estern part of the horst structure. All the large horizontal 
faults in the ac-direction ( ’’antiapenninic” ), seen in the 
aerial photos, are more difficult to recognize on the 
LANDSAT image , 

The existence of two different b-axes ( see chapter : "Examples 
'of Characteristic Pole Diagrams out of the Test Site") indi- 
ca'tes two deformation stages, the first of which, at least, is 
active still in the upper Trias sic, but no longer is active 

I 

in the Oligocene. A multistage deformation is postulated by 
several scientists investigating the Geology of the Apennines 
j(Baldacci, et al, , 1967, Bortolotti, 1966, Giannini, 1962, 
Sestini, 1970 , Reutter, 1968). A possible explanation for the 
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' difference of the two b-directions is given by Boccaletti, 
Giiazzone ( 1968, 1970 ). They postulate an anti-clockwise 

‘ I 1 

rotation of the Italian Peninsula, which brought the axes 
to their present directions. The structxaral data reported 
by these authors ,( e.g. the E-V/ stress direction in the 
Tuscany, corresponding to N-S fold axes ) cannot be con- 
firmed in oxir mapping area, Wunderlich’s (1966) opinion 
that the stress direction rotated counter-clockwise to the 
■ WKVf-ESE trend on the margin of the Po plain is better 
suited| to the structural facts. A new concept was proposed 
' by • our- research team: based on the interpretation of LAJIDSAT- 
j- 1 images we contend that a clockwise rotation of a block 
formed Tuscany and part of Umbria. Such a rotation results 
from a lefthanded horizontal fault following the Arno valley 
.•■and another N-S striking left handed horizontal fault between 
<the towns Ancona and Latina ( see previous part of this 
.report). ‘ 

J'i 

5 - 3 


i ’ I 

The relations between field work and aerial photointerpreta- 
■ tion are demonstrated in figureig and 2o. Figure 19 shows the 
main lineaments touching the Mesozoic ridge in the area 
around Trequanda, as seen on aerial photographs. 

Recognizable displacements are marked as well. Figvire 20 is a 
structure map of the study area S of Trequanda. This map was 
prepared from data supplied by intensive field Investigations 
which paid special attention to large lineaments visible bn 
aerial and satellite images. Such combined work applied to a 
small ground area has provided a solid base for the inter-re.- 
gional interpretation of jLANDSAT images ( as performed in the 
previous chapter). The combination of field measurement and 
lineament analysis indicated the existence of two. fabric 
systems with different orientation and suggested the concept 
of cloclm-/ise rotation of Tuscany-Umbria. ' 



Fig. 21: Rose diagrams 
of the phot olineat ions 
of the test aerea 




Pliozjin 


242,2 km 


Gesamtgchi ot 


?.rL.,4 km 



ZENTRALSTEULE FUR GEO - PHOTOQRAMMETRIE UNO FERNERKUNDUNG 


8LATT 32 


The Photo Lineaments 

The lineament map (Figure 19 ) has been statistically 
evaluated in three separate zones : 

zone 1 : formed of Pliocene sands 
zone 2; formed of Mesozoic limestones and cherty- 
calcareous-marly alternating strata 
zone 3i formed of Macigno sandstones and Argii'iHe 
Scagliose 

The rose~diagrains show the percentage of all the lineaments 
for 10*^“Classes. Looking at the four rose-diagrams (Fig. 21) (one 
for each specific zone and one for the whole region ) the 
most striking feature is the similarity . of the diagi-’ams 
for the whole region and for the Pliocene, This is caused 
by the fact that Pliocene in this region covers by far 
the largest areas. So it has great influence on the percen- 
tage of lineaments of the whole region. The evaluation of 
the three zones agrees with the conclusions reached by the 
interpretation of equal area projection diagrams of joint 
poles. 

In the Mesozoic limestones the early deformation stage is 
reflected by lineaments in ac,^ -direction accompanied by 
hko^ shear directions. The stratification was nearly hori- 
zontal during this early deformation stage. The yotmger 
stage (during upfolding) is represented by aC 2 and bC 2 . 

Very important are E-W-striking hk 02 -lineaments reaching 
the absolute maximum in this zone. Zone 2 (Agrille Scagliose 
and Macigno) show especially. ac 2 > bc 2 and hko 2 formed by 
the younger deformation stage. The scattering of the maxima 
can be explained by the different competences of arenaceous 
and argillaceous beds. A last stage with 125° striking be- 
traces seems to have been important in late orogenic 
block -faulting. In zone ? (transgressive Pliocene sediments) 
the fracture traces are not statistically distributed, but 
the directions of the undei lying Mesozoic rocks seem to 
penetrate through the youngei- sedimentary covering. ■ 
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5.4.Some_Methodologlcal_As2ects_on_the__Evaluation_of ■ ' 
s o “Te c t onical _M§§sur ement s 

' jAproximately 15oo joints and 5oo bedding planes were 
jmeasTired in the field. These measurements were plotted 
in pole~diagrams (equal area projections), Interpre- 

■ tat ion of the diagrams is based on theoretic.al principles 
developed by; R. Adler et al, through intensive -research 
in the' -coal, district of the Rtohr Region, Western Germany. 
Adler found, by the measurement of several hundred points, 
that two deformation systems, corresponding, to two defor- 
mation stages, always .existed. In the younger -stage, the 
older joints and faults are only revived, forming suitable 
anglers to the new main stress direction. In the other case, 

.totally new parting planes have originated. In addition to 
*ijhese two systems in most cases a third system is recogniz- 
abie, caused by a special stress field in the s-urro\mdings 
of .large scale faults. One aspect of the field work was' 

^ 'to verify the influence o'f such large lineaments on the 

I 

tectonic deformation system, as can be seen in- LAiJDSAT- 
Images. 

I ... ... 

Figure 22 shows an ideal stress-cube with the correspond- 
ing positions of the dividing surfaces. Figure 23 illustrates 
^ surface appears as a great circle in the equal area 
, projection diagram (Schmidt 'Wet) . For statistical evalu- 
..'ktijons the use- of the projection poles of the surfaces is 

■ more suitable since numerous great circles on a single 
[diagram are disturbing to the interpreter. 

The following symbols are used in the pole-diagrams; 

•• .1 * ! 


o 

ac 


.A « hko 

o 

=B be 


o A= hoi 

0 

» ab 


O A=Ohl 


" ’ J 
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The following pole-diagrams ( Fig. 24 - 26 ) are only 3 
characteristic examples for the tectonic data, which were 
taken directly or photogrammatrically. 


The ooint-pole diagrams ( Fig. 24 , 25)clearly show 

the -existence of two deformation systems, characterized 
by typical configuration of parting planes, as follows: 

System 1 is characterized by well -developed ac-fractxires, 
whereas be and hko are less important. Furthermore, hoi 
and okl can be seen. The b-axes of the various diagrams 
have a strike between 166*^ and 17o°. This fabric can 
surely be attributed to the stratification (ab great-circle 
representing bedding planes) and for this reason it must 
have been formed before folding. 

On the contrary, the second fabric system is composed 
mainly of ab- jointing. Ac and hko are represented less 
or are totally lacking. On the other hand, hoi and ab 
(in the massive Liassic limestones) are well-developed. 

The b 2 ~axes sho\if strike directions of 142-158*^. They are 
identical with the fi -axes of the bedding -pole diagrams, ( ?ig,g» 3 ) 
The aDove-mentioned characteristics of the b 2 -axes and ■ 
the significant assemb3,ing of joints Indicate the origin 
of the second fabric in a later deformation stage, during 
upholding of the strata. 

In addition to those two basic deformation systems, a 
special fabric can be related to important fault structures 
(thrust faults and horizontal faults). Such systems are 
defined b^)’ secondary shear jointing, as described in the 
introductory chapter. 
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The evaluation of the tectonics of the Northern and. 
Central Apennines as described in the previous 
chapters is based on field work, interpretation of 
aerial photographs and satellite imagery. Aerial 
photographs and especially field work cover small 
areas. Therefore only, the synoptic view of the- 
satellite images allowed to gain the results, which 
partly lead to new ideas on the tectonic systems of 
the Apennines. 

The methods applied dviring this work are not only of 
scientific value but also give a basis for the ex- 
ploration of tectonic controlled mineral deposits by 
satellite imagery. 
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B) GEO-TECTONIC STUDIES IN EASTERN ALPS 
INTRODUCTION AND METHOD 

A section across the eastern Alps has been analysed 
using the imagery from the LANDSAT-1 and 2. It has been 
shown, as a result of the present studies, that's great 
potential of these satellite images lies in delineating 
active movement zones of today. Presence of such move- 
ment zones along with their sense of movement can be 
deciphered, on the satellite images, by very aceuxats 
and detailed mapping of lineations and study of deformation 
of landforra structures and river courses. 

Conventional methods of photo-interpretation utilizing 
texture, tone and morphology, have been applied to these imapes 
for delineating faults and lineamentsj^:‘^Special attention wac 
paid to the trends and their deformation, if any, of the 3.atnijc.vm 
structures and river courses which has helped In deciphering 
direction of movements along the faults.' For examples; study 
of the courses of tributaries of the Adige river in the Bolzano- 
Kerano region has indicated presence of an active fault (' sinis- 
tral sense of movement ) along the river valley. The wain Adige 
river flows NW-SE in this tract. Tributaries on the eastern side of 
the Adige river turn southwards while approaching the main river, 
whereas those on the western side turn northwards, thus simula- 
ting. a typical drag-effect ( Fig.29).Tl^e southwards turn of the 
tributaries on the eastern side may be accounted for also due 
to general loss of elevation southwards. However, the fact that 
the tributaries on the western side turn anomalously northwards — 
towards the region of general higher elevation and against the 
direction of flow of the main river - brings clearly out the 
presence of active movements along the Adige river valley in 


This work has been carried out by Dr.R.P Gupta, guesi 
scientist from the Department of Geology and Geophysics, 
University of Roorkee, India. 
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Ficf -27; Liaeation tnap of the eastern Alps as inter- 
preted from the LAKDSAT-"* and 2 ima°-es 
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Fig . 29:TurninE; of che tributaries indicnring ■'c iv- si :;st* i 

faulting rlong the Adige river valley in ’hi: B)- 

Merano CM) region. 

LANDSAT-1 , dated channel 7 

this area, the sense of movement being sinistral. Tne lineament 
can be traced and extended in the mountainous tract on either 
Bide, on the imagery, where also movements with sinistral 
sense along the lineament are confirmed. It may be added here 
that as the satellite images provide vertical views, only hori- 
zontal movements can be deciphered on the images. However, the 
swampy areas occurring in the south-western part of the river 
valley in consideration indicate that the possibility of con- 
currant vertical movements with western block as downthrow side, 
is also not ruled out. Similarly, accurate detail mapping 
of phot o-lineations representing geomorphological features 
has shown the presence of drag-ef f e .‘ts pertaining to dextral 
movements along the Insubric-Tonale Li ne ( Fig. 3q) . 

REPRODUCIBILITY OP THli 
ORIGINAL PAGE IS POOR 
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Fipr. 30: Drag effects a lung the Insubrj c-Ton xle Line suggesting 
movements with dextral sense. 

LANDS * T -1 , dated 7.1o.l972, “ 0,P • ■* , ’ 
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Interpretations have been made by superimp^osing observa- 
tions on raulti—time imagery. Images obtained on 13,8,72, 31.8,72 
and 7.10.72 by LANDSAT-1 and those on 7.8.75 by LANDSAT-2 ' 

have been made use of. This, on the one hand, has nelped ^n 
removing errors due to illumination conditions, .:such as 
position and angle of sun leading to shadows. •( Bodechtel 

and Nithack, 1974 ) , and on the other hand, has strengrhened 
validity of the observations and interpretations therefrom. 


RESULTS AND DISCUSSION 

It is apparent ( Fig.2J) that nearly all lineaments are 

characterized by several long to small en-echelon minor linea- 

tions constituting a zone. So, basically we have to visualize 

zones rather than planes such that each zone is comprised of 

numerous small to large paraYlely arranged cn-echelon minor zones 

snd 'lurfaces. Fig. 28 shows a generalized map with major tectonic 

zones and lineaments. A few of the lineaments are seer, to be 

very extensive, being traceble from one end to other end of the 

investigated area and running across the entire section of the 

Alps. Some of the lineaments have been off-set by others. Besides, 

as mentioned earlier, several lineaments have been found to be 

active(Fig. ' 39 .,So,?i:?)Results obtained in the present investiga tioiis 

render it necessary to make a special mention,, among many others, 

c 

of the Insubric-Tonale Line and Pustertal Line — the two rscgnized 
major geotectonic boundaries in the eastern Alps ( Fig.2j^ 

Insubric-Tonale Line :0n the satellite images, this lineament 
appears to consist of a number of minor overlapping iineations 
and shows displacements along other intersecting presumably later 
lineaments (Fig..27 ,-2£)The present investigations, as shown (Fig.^^o, 
indicate m movements (dextral) along this Line. 

It. is worthwhile to recall that Gansser(1968 ,p. 139) on the basis 
of field investigations in this region wrote about temptations 
to postulate lateral displacements and Laubschsr (1971) based on 
palinspastic restorations involving the Dinaride Alps and the 
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Fig. 31: Regional tectonic map 
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Apennines inferred 3oo km. of post-Oligocene slip, both suggesting 

similar sense of displacement as deduced above ( also c.f. Laubscher, 
a ,b 

1973; Ernst ,1975^ .Further, the Insubric-Tonale Line is commonly shown 
to be abruptly terminated on the east by the Giudicaria Line (Fig .31) 

However, the present investigations indicate continuation, of this 
trend in the acidic extrusives lying on the eastern sidie of the 
Giudicaria Line and the Adige river(Fig. 27) Though the .intensity 
of development of this lineament in the acidic extrusives 
is not so pronounced, the mere presence is striking. In the 
geological past, the Insubric-Tonale Line, of course, seeras 
to have been cut-short (or displaced) by the Giudicaria Lins. 
However, later movements along the Insubric-Tonale must be 
postulated to explain its eastward continuation in the acidic 
extrusives. Further, though the seismic investigations have 
failed to record movements along this Line ( in Gansser, 1968)', 
the results of xhe present investigations that all along this 
Line the georaorphological features are aligned simulating drag 
effects, strongly suggest that this geotectonic boundary is still 
aetive or at least has been active in the recent past. 

Pus'tdrtal Line ; Similar to the Insubric-Tonale Line, the Pustertal 
Line is a major geotectonic feature(Fig. 31) The geological 
field mapping shows no indication of continuation of this Line 
into the area of present study and therefore it has been, by and 
large, considered to bs abruptly terminated by the Giudicaria Line 
on the west. Continuation of this Line westwards across the 
Giudicaria Line was first detected on the ERTS-1 imagery by 
Bodechtel and Lajnmsr8r(1973) . In the present study it has been 
further established that this Line constitutes a zone of active 
movements as accurate mapping of the georaorphological features 
adjoining this lineament on the imagery • . . 

has revealed drag-effects with dextral sense of movement. Its 
westward continuation across the Giudicaria Lins and nresence of 

• V 

drag-effects indicate its later rejuvenation and still active 
nature (dextral). It may be of interest here to mention that 
on the basis of 'repetitive precise levelling measurements over 
time, active vertical movements along this zone at the rate of 
nearly 1 ram. /year have been found by Senftle and Exner (19'73) . 
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system of Post-Alpine Faults ~ Relation to the Present Day 
Stress Field 


An interesting feature is the presence of several VI-E to 
NW-SE trending lineaments some of which have "been neither 
mapped so far, nor could he visualized earlier (Fig, 

On the satellite imagery, they appear as rather unpromine'n,t tine 
lineaments which failed to make an impression during the. first 
hand study for major lineaments. It was only during detailed 
analysis that their presence became clear and thenceforth has 
been confirmed on separate observations on both LANDSAT-1 and 2 
images. It is important to note that it is not the case of an 
isolatory occurrence; on the contrary they appear to form a set 
and hence certainly demand attention .These lineaments are fine 
but extensive and cut across all earlier trends and boundaries 
thus indicating that they are post-Alpine in age. Often along 
their length they are marked with geomorphological evidences, 
in the form of drag-effects, indicating movements of evidently 
recent nature. Some of these lineaments partly follow the old 
major tectonic zones, like the Insubric-Tonale Line and the 
Pustertal Line, apparently extending and rejuvenating them. 
Further, a relation between these lineaments and the stress field 
leading to folding and uplift of the Alps (maximum principal 
stress at K 15-20^ in this area , see page 51-53} is difficult to 
arrive at. On the other hand, it is easy no link these lineaments 
to the . stress field of the present times. On the basis 

of in-situ stress measurements and fault plane studies of earth- 
quakes, maximum principal stress of the present day stress field 
in central .Europe has been found to be "W-E to WNW-ESE (Ritsema , 
1969; lilies, 1974,1975; Ranalli and Chandler, 1975), which 
taking the intermediate principal stress to be vertical, would 
fairly well explain the development of these post— Alpine seemingly 
active zones of^ movements (Fig. 32) . 


WBODUCIBlIiTy 
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Sect ion studied . 


Post- A (pine movement zones with ^ 
sense of horizontal displacement 


Direction of the maximum' principal 
stress of the present day stress field 
jnCentralEufope 

Fig.t 32: Relation between post-Alpine movement zones 
and the present day stress field 
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Predominant Llneation Trends - Relation to Stress Field 


A trend diagram depicting relative predominance of various 
lineation trends observed on, the satellite images of the area, is 

3 T • , • ’ 

shown in Fig._^3. Clearly three major directions sitand out:(i) 

N 45-225, (S^); (ii) N 15-195, (T ); and (iii) K 345-165, (S^) . 

Out of these three lineation groups, the two S- and 

O X 2 * 

oriented at N 45 and N 345, respectively are considered 

to be of shear origin and the third, T, trending N 15, to 

be of tensile origin (all the three being of oogenetic 

nature) on the basis of following arguraentst 

(i) A number of movements planes are observed , orienxed 

£ 

at nearly N 45 (fig. 27 and hence the- lineation group 
is of shear origin. 

(ii) The lineation groups and appear to form a conjugate ~ 

set. Here a few additional remarks seem warranted. The pile of 

rocks involved here range from rather incompetent sedimentary 

deposits with surficial loose moraines to low to high grade meta- 

morphics and crystallines(Fig. '??,) .Regional geology after Carte 

tectonique Internationale de I'Europe, l:2,5oo,ooo) . Evidently, 

they have wide diversity in their failure characteristics. 

♦ 

However, a corollary to the Coulomb-Navi er criterion for 
shear failure widely applied to rocks (c.f. Anderson ,1951; 

Farmer, 1968; Price ,* 1959 ,1966) states that under normal 
stress conditions the dihedral angle of the shear fractures 
( 26) is related to the coefficient of internal friction ^ 
as follows; 




Ian 

e * 
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Fig. ?3: Trand diagram of Liueations 
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• Fig, 34: Preferential sinistral movements 
on the NE-SW trending lineaments 

'i 


ESPROBUCIBlLrrY OF THE 
OBIGMAL ?AGB IS POOR- 
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For most natural rock substances the angle of coefficient 
of internal friction (ji) being around 3o— 45 , the dihedral 
angle ( 26 ) ranges around 45-60 , 'which is in conformity 
to the angle between and Sg C taking intermediate 
principal stress to be vertical ) . Further, the above 
relation bet'ween the dihedral angle and the coefficient 
of internal friction implies that the dihedral angle 
increases with decreasing frictional resistance. This fea- 
ture seems to be corroborated in the present study as the 
Sx and S2 lineations seem to widen out towards the northern 
tract from the central parts of the Alps (fig, g.7) .This 
’’refraction” observed in this area, is attributable to the - 

decreasing competency of the medium, as the northern tract 

C 

is underlain by largely molasse and surfi^l moraine deposits 
and the central mountainous range is underlain by coyipexent 
sedimentary and crystallines rocks. 

(iii) A third set of lineament (T) is ori.ented at N 15. Few 
displacements are noted parallel to this set. On the other 
hand, many of the larger lakes are observed to lie along this 
set of lineament, for example, Starnberger Lake and Ammer 

Lake in the North and Garda Lake (?) and Como Lake in the South 

t c 

(fig. 27) .Further , this trend (N 15 ) coincides with the di- 
rection of maximum principal stress as derived from the 
orientation of and Sg conjugate .shear sets . Therefore, 
it seems logical to conclude that the lineament set T has 
cogenetically developed as a result of tensile failure. 

Hence, the above interpretations that the S-, S2 T 

groups of lineaments have, cogenstically developed and that 

and S2 are of shear origin and T of tensile origin, are well- 
interknit and appear to be logical. It is interesting to find 
how the direction of ma-xrmun: orinclpal sty.ass deduced 

as above from the lineament analysis of Alpine section 
compares with that independently inferred for the adjoi- 
ning areas by ether workers. A major ' geotectonic 
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feature adjoining the Alps is the Rhinegrahen system, con- 
sidered to be more or less contemporaneous with the Alpine 
orogeny. The Rhinegrahen trending NNE-SSV/, evidences crustal 
dilation normal to the graben axis, of about 4.8 km. and an 

average direction of maximum principal stress of abogfii N 2o ( INHE) 
must be concluded for the beginning of rifring in Mi-tfu'le Eocene 
times (lilies, 1974, 1975 ). This is in striking conformity to 
the direction deduced as N 15 from the above lineament analysis 
f®^ the Alpine section .Further , the above inferences regarding average 
direction of maximum principal stress in late Mesozoic - early 
Tertiary times from the macrotectonic observations are corroborated 
by the results of statistical evaluation of field measurements 
on horizontal stylolites in unfolded Mesozoic platform setilments , 
made all over Central Europe ( Beiersdorf, Wagner, Wunderlich 
and Plessmann,in lilies, op.cit. ) 


The presently concluded situation that the then prevalent 
maximum principal stress seems to have been not perpendicular 
(and hence the minimum principal stress not parallel) to the ■ 
general trend of the rocks in this Alpine section needs 
explanation.lt seems that this abnormality^ been caused by the 
presence of a hindering block on the south-east. Th^lternative 
possibility is by postulating two different phases of deformation 
the first phase (with average maximum principal stress trending 
N 35o) bringing about the regional folding and uplift of the 
Alps, followed by a second phase (with average maximum principal 
stress oriented at N 15-20) causing largely translational 
movements, shearing and faulting but retaining the pre-existing 
strike trends. The author tends to favour the first possibility 
as separate phases of folding and faulting are difficult to 
visualize. Often the two phenomena' are found to be^ concurran-t and 
closely associated with each-other. If the main phase 
of deformation had maximum principal stress at N 350 , the 


related tectonic impressions which would have been of the Same 
degree of magnitude and intensity, should have been visible 
and thus borne evidence to that. Further, uplift of the Alps 
associated with northwards movement of the Adriatic block 
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has been long postulated ( foi* example , Aubouin , 1964) and the 
same block, which lies to the south-east of the area, could 
have also provided the hindering resistance referred to 
above. Besides, the presence of such a block would be further 
in conformity ' to the observed predominance of NE-SW (S^) 
lineations with sinistral sense of horizontal displac.ament , 
described below. \ 


Lineations and the Giudicaria Line 


A striking feature is the large number of NE-SW trending 
lineations ( S^) with predominantly sinistral sense of movement, 
observed throughout the region of study (Fig. 27) .They are 
developed with varying intensity in the area .and are more 
prominently developed in the vicinity of Giudicaria Line. 

From field investigations Trevisan ( in Laubscher^ 1972) 
reported the large system of NE-SW trending sinistral faults 
in Giudicaria. Evidently , intense development of set of 
shear planes and zones seems to have controlled the development 
of Giudicaria Line, 

The reasons why movements have taken place « largely on the S^-set 
of shear surfaces and not on the Sg-set and that too with 
predominantly sinistral sense are also" clear. The general strike 

O 

of rocks is N 75-8o (Fig. 27) .The maximum principal stress has , 

* £. 

been concluded to be directed at N 15, giving the conjugate set 

t ^ 

of shear surfaces at N 45 (S^) and N 345 (Sg) - Thus, the shear 

surfaces S„ trend nearly perpendicular t6 the general strike 
^ ■ c 

of rocks whereas those of S^-set make acute angle (30-35) with 
the latter(Fig.S^) .In such a state, movements on the S.j^- group 
of surfaces ,'withjsinistral sense, would be evidently preferred. 
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CONCLUSIONS 

On the basis of above studies on LANDSAT--1 and 2 images of 
eastern Alps, the following broad conclusions may be drawn : 

(a) A number of W-E to NW-SE trending post-Alpine active 
lineaments are present. These lineaments can be logically 
related to the present-day stress field. 

(b) There are three dominant lineation-sets in the eastern 
Alps. They have co-genetically developed as a result of 
shear and tensile failures, due to stresses with maximum 

c 

principal stress averagely trending N 15. 

(c) The Giudicaria Line seems to have been controlled by the 
NE-SW trending lineaments over which -movements with sinistral 
sense have been preferred. 

It goes without saying that as the satellite images provide very 
accurate map of ground features, a great potential of these 
images lies in their detailed and accurate mapping. This can 
unearth valuable information of local and regional importance, 
unobtainable from any other source. 
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C APPLICATION OF LANDS AT DATA FOR LAND USE MAPPING BY 
CONVENTIONAL INTERPRETATION -TECHNIQUES 


For many earth scientific activities,' land use maps Sorm 
the basis for further planning and inventory prepara'tdLon. 
Very often the existing land use maps are more than 2c to 
4o years old and do not correspond any longer to the 
existing situation. 

In order to prepare a land use map on a conventional basis, 
many man-years are necessary to conduct intensive field 
work and supporting aerial photo interpretation. The Italian 
authorities are faced v/ith the problem of producing nev/ 
land use maps of the mountainous regions on a scale of 
1 : 5oo’ooo. Especially in mountainous regions, the cost 
factor has always been a severe constraint, hampering 
completion of this project. Preliminary evaluations of 
LANDSAT-data, has indicated the high potential of multi- 
spectral satellite data and its suitability for land use 
mapping. On the of our background experience in 

aerial photointerpretation and satellite data evaluation 
technia_ues, v/e approached this problem by determining the 
feasibility of vegetation mapping through application of 
conventional interpretation techniques to LANDSAT-data. The 
tasks seem to be solvable with respect to the small scale 
and the few main or level one categori'e-S that have to be 
determined. The different categories of interest are: 

1. barren land (rocks, gravel) 

2. permanent snow cover, glaciers 

3. tarban areas (including industrial areas and larger 
■ transportation patt^-rns - streets) 

4- pastureland (areas covered by grass) 

5. agricultural areas (including plantations) 
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6. forested areas 

7. water bodies, rivers and lakes 


1 1 Data Source 

LAIJDSAT-1 MSS images taken, during ‘72 and *73 wei^iiUssd 
as the data source* The nine inch negatives (scal%f1: 
i.ooo.ooo) were photographically enlarged on a scala ot 
1:25o.ooo. For most of the areas the information loss . 
acquired through repeated photographic processing of the 
images did not interfere v^ith the detection of level 1 
categories. For some key-areas a direct production of 
images from tape were carried out by PR AKLA- Company. 

These images, showing maximum information of the MSS 
data, were generated on a scale of 1:2oo.ooo. Furthermore 
for some middle Italian and Sicilian areas Skylab photo- 
graphic data in panchromatic and false colour were pro- 
duced. By this means a total cloud-free-coverage of the 
areas of interest was obtained. 


2. Evaluation Techniques 

The basis of any photointerpretation is not onl 3 r the 
experience of the interpreter, but also detailed background 
information. The background information for our project 
consisted of; 

a) existing land use maps on a scale of 1:2oo.ooo, 
produced 2o - 4o years before 

b) existing field data obtained through intensive 
field Investigations of the last ten years 

c) detailed vegetation profiles. 

Information gained through field work proved to be the 
most actual information and could be utilized with high 
confidence. 


RTOODUcmiLiTr Of fis 

ORIGINAL PAGE B 
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2 . 

2 .. 2,1 LANDSAT-1 multi spectral scanner data 

For conventional photointerpretation there are no sig- 
nificant differences of inf ormation between MSS 4 and 5 
and between MSS 6 and 7. 

On the other hand MSS 5 and 7 sho\ir the best contrast of 
grey tones in vegetated and non-vegetated areas-. There- 
fore images of these two bands were chosen for the land- 
use mapping. 

From an earth scientific point of viev/', the utilization of 
data representing different vegetation cycles would be the 
optimum. 

This requirement could not be ftilfilled over all areas 
because of insufficient continuous data. It turned out that for 
the production of a preliminary base map, channel 5 was op- 
timally suited. Band 7 was utilized in order to obtain 
additional information for discriminating various types 
of surface features. So, for e 3 cample , the use of band 7 
allows for the separation of coniferous forest from mixed 
forest, the better delineation of agricultural areas and 
the clear identification of barren land. For the level 1 
classification the sensitivity of band 7 to soil moisture 
has been somev/hat disadvantageous. 

The diagrams of figure35 represent a qualitative comparison 
of grey tone classes obtained from band 7 and 5. In the co- *» 
ordinate-system the land use' categories were plotted as a 
function of their grey tone in both channels. It is evident, 
that such schematic representations only gives a very rough 
idea about the principles of evaluation. Of course such 
diagrams do not reflect the experience -of .the interpreter, 
his ability to fifferentiate subtle grey tones or his 
analytic ability under various environmental conditions. 
Furthermore both diagrams do not express important criteria 
such as texture and morphology and also the synoptic know- 
ledge of the geological conditions. An interpretative -- 
decision can be made by a combined analysis of all criteria, 
as shown for grass in figtire 33 . 
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I dentl f Ic.at ion of grass lane] 

Schematic diagram, for the I nterpreta t ion of L^MDSAT WSS-data 
(Legend: GT = grey tone, L v ligKt, LL « very Ught, 

M = medium, D‘ = dark) 
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2*2,2 Skylab data (colour and false colour images) 

Colour photographs of Skylah (for example 19o a, taken in 
Septemher 1975), were projected on- a scale of 25o.ooo. Due 
to the outstanding colour quality of these images, the 
following surface features could be clearly identified; 
forested areas, agriculturally used areas, grassland and 
barren land. False coloxir coir.posits of Skylab dataand also • 
of LAKDSAT-1 MSS data served as basic information for ■ - 
differentiation areas with grassland from agricultural 
use regions. The determination of borders between those 
two features has been very difficult with black and white 
images taken during late autumn, summer or winter. Curing 
these seasons, the fields are not planted or are already 
harvested. In relatively dry areas grassland and agri- 
cultural areas appear in relatively vuiidentifiable grey 
tones •v^?hen the field patterns are far below the resolution 
capabilities of the sensor system. The only reliable. cri- 
terion for differentiation is a distinct red colour re- 
presented in false colour images. 

2.5 Gra^^ Tone and Colour Tone 

In respect to the physical properties, various surface 
features reflect incoming radiation in different inten- 
sity levels. The reflection Intensity is translated on 
photographic data into a gray or coloured tone. Therefore 

the -Dhotographic density is one of the most important inter 
pretative criterion. In addition, the density is very 

sensitive to uncontrollable factors, caused by natural 
phenomena or by photographic processing techniqtxes. 
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2.3.1 . The_l22act_of_^atural_Parameters_on_'the_Photo^a2hic 
Density 

In general pictorial representations of our earth surface 
are a comhination of all possible factors of v/hich the 
image density only pla 3 Fs one part. Problems aris^ycon-* 
cerning the identification or surface features, regarded 
to be homogeneous from any earth scientific point of vie^;. 
This is due to the fact that such a homogenity does not 
necessaril 3 '‘ correspond to a homogeneous distribution of 
gray tones. Therefore in order to decide whether to deal 
v/ith category a or b, further decision-malcing criteria 
have to be introduced. So for example, a well cultivated 
fertilized medixim humid meadow may have a darker grey tone 
than an imcultivated grassland under the same environmental 
conditions. Grassland in an Alpine valley, for example, 
has a different grey value from grassland in a south 
Italian Apennine valley. Grass growing on clay has a 
different dark grey tone than grass on limestone. The grey 
tone, the intensity of reflections, is mainljF dependent 
on natural factors such as soil moisture, soil type, cli- 
matic zone and local climate. The image, figure '37, re-* 
presents some of the basic relationships betv;een environ- 
mental parameters and grey values. Another important 
factor, not taken into account, is the exposure of the 
surface to an illumination source - the Sun. In areas 
with great relief the identification of surface features 
is extremely difficult, not only due to the shadows 

themselves, but elso to the change of the relative sunangle. 
Another Important factor is the changing state of vege- 
tation during the year, i.e. the change of reflectance 
depending on the seasons. Therefore multitime imagery 
allows better identification and further differentiation 
of land-use categories, 

2.3.2. Impact_of_ Tecnical_Fac^^ PhotojrsroM 

The main technical factors influencing the density of 
images are related to the image generation itself. 
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Fig. 37 Grey-tone levels of grass 

Schematic representation for the factors soil, 
soil-moisture and climatic zone 

Concerning the generation of images from a digital tape, 
restrictions result from the dynamic range of the output 
devise. The equipment used for data generation has an 
overall dynamic range of 16 grey levels. In comparison 
vrith the theoretical p erf or mane e . of the MSS scanner, with 
a dynamic range of 128 (MSS 5> and 64 (MSS 7.) grey levels, 
the applied generated images represent reduced informa- 
tion. In order to minimize this effect the available inter- 
val of 16 grey tones was adapted to the actual sensitivity 
range of the scenes. In order to obtain maximtun grey reso- 
lution for the areas of interest-, 'the adaption was' carried 
out on scenes that did not show maximum reflective contrasts, 
such as snow and watei-^ in infrared images. A second important 
parameter which influences the quality of the film material 
is the degradation effect resulting from photolab processings 
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FiEr.3 8 

Identification of wooded areas 

Schenatic diagram for [jhe interpretation of LANOSAT MSS* data 

(Legend; GT ¥ grey tone,! L •= light, LL - very light, 

M = medium, D i| dark, DD = very dark) 


BliPRODXJCIBrLlTY OF THE 
S^n-TNA-L^PAaB IS. PQOB-^. 
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such as enlarging etc. During these processes the grey 
tones are influenced by the various film emulsions f by 
various exposxire times and by the state of the chemical 
ingredients. 

Another degradation effect which occurs mainly is chamiel 
7 is caused by the high contrast between v/ater end llsaa 
masses. The extreme dark water areas interfere in . s ' 
v;ay that the coastal areas always appear in darken- • 

tones. The above-mentioned photographic factors may es- 
pecially degrade the possibilities and the accuracy of 
conventional photointerpretation techniques. In other words 
the interpreter has to revise his interpretation results 
not only from image to image but also v/ithin the evaluation 
of one scene. 

3, Possibilities and Limitations of the Conventional 
Classification 

The possibilities of identifying surface features is 
extremely dependant on the resolution capabilities of the 
sensor and on the target size. Figure 38 demonstrates to 
what degree under optimum conditions, the .differentiation 
of trees can be made. Vfith. regard to trees or tree-cones, 
two points have to be taken into consideration. l*/ith 
respect to the resolving power to the multispectral 
seamier, larger areas v/ith a dense tree population can 
be identified as tmiform reflecting areas. It is under- 
stood that with a dense population of trees there is no 
interfering reflection from the ground. Furthermore be- 
* cause of the poor resolution of the satellite data, the 
differentiation betvFeeu bushes and forests, can not be 
made. 

Furthermore conventional definition based on the percen- 
tage of trees within a certain area can not be made. 

To a certain degree, grey tone and background Information 
enable the interpreter to map areas with different densi- 
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ties of wood. But the geometric resolution of LA^TBSAT-"^ 
images does not allow to map woods by quantitative measure- 
ments of percentage of trees per area. Therefore areas 
with for example lo % trees, which are conventionally de- 
fined as wood do not appear on the evaluated land-use map. 

Further difficulties arise when comparing forsts and '• 
plantations. Densely planted hazels agrumenplantages 
and vineyards have the same reflection characteristic ’k’s 
deciduous forest, at least for the eye, A basic local 
knowledge and the experience in examing textural infor- 
mation may sometimes help to overcome above-mentioned 
difficulties. 

Plantations, for example, are very often characterized 
by geometric borders, and in contrast, vineyards often 
show a kind of chaotic texture. 

Cultivations with clearly separated rows of plants such 
as oliveyards are characterized by a uniform medium grey 
value. These grey values are originated by intensities 
reflected from the ground coverage and the foliage of 
the trees. The ' following images demonstrate areas with 
plant units which have characteristic signatures. Figure 39 
shows coniferous and deciduous forest near Pratomagno, 
mountains south east of Florence, figure 4o mediterranean 
oaks in the Apennine west of Rieti , figure 41 oliveyards 
of the Tuscany, figure 42 an area with bush and grassland 
near Valle Taleggio (Begamasc Alps) . 

Especially the last example (fig. 39) indicates the diffi- 
culties which result from the poor geometric and spectral 
resolution of LANDSAT-l-images. 

Areas with mixed reflections are very difficult to classify 
without knov/ledge of the local condition.^. Nor-vnally in such 
cases textures can not be detected any more and the surfaces 
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appear rather uniform. The limitations of differentiating 
between various surface features becomes evident if we 
are looking at interface regions betv/een meadovrs and 
areas with dried-out grass and barren rocks. Barren rocks 
are clearly detectable but a distinct tracing of the 
border between the grass land and areas without vegetation 
is impossible. 
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Pratonagno Vountains (detail from LANDSAT-1 image of 12. Aug. 1972) 

a) MSS 5 (dark - wooded areas) 

b) M*^S 7; light = deciduous wood ( L - in general, K - chestnut ) 

d^rk - coniferous ( N ) 

c) Ground truth photograph, showing coniferous (dark) surrounded 
by deciduous wood 
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Fig:. 40 : 

Mountain ranges west of Rieti (detail from LAND'^AT-l image 
of 6. Feb. 1972) 

a) MSS 5 (dark - mediterranean oaks) 

b) Ground truth photograph taken in June 1973 near Amelia. 
Dark - mediterranean oaks 
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Fig. 41 ; j 

Lower valley of the Corrla 
River, east of Piowblao 
(detail from LAFDSAT-1 image 
of 12 ,‘ug. 1972) 

MSS 5; dark - woods 

middle grey - olive yard 
(arrow) 


Fig. 42 ; 

Valle Taleggio (Bergamasc 
Alps) 

LANDSAT-1, MSS 5, 7. Oct.? 972 

(a) does not show the distinct 
vegetation pattern as the 
ground truth photograph (b) 
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Discussion of the Results 


The possihllities of conventional interpretation mentioned 
above base on extreme situations. With respect to the 7 
level I categories the results are good. 

The land-use classification followed the USGS- system which 
ANDERSON, HARDY & ROACH (1972)* established for the inter- 
pretation of remote sensing data. 

For our special problems there was introduced the level I 
category "gras^and”, which in the USGS-system belongs to 
•'agricultural land", if it is pasture, and to "range land" 
if it is uncultivated. ; 

Th.e following table (Fig. 43) -shows to which degree LANDSAT 
and SKYLAB images can be used for mapping land-use categories 
of level I and II, The table is based on our' experiences 
during this work and therefore is valuable for European 
moiantainous regions. 

j Good satellite dataand good background information enable 
the interprete to subdived the level I categories no. 2, 4, 

5 ,^d 6 in level III and IV categories. The following 
schemes show these subdivisions, which are not included in 
the above USGS-system, ^ 


* J.R. ANDERSON, E.E. HARDY & J. T. ROACH: A Land-Use Classi- 
fication System for Use with Remote Sensing Data. - Geol, 
Surv. Circular, 671 . 1-16, Washington 1972 
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Fig, : Quality of Landsat- and Skylab-iraagery for land-use mapping if 
level I and II categories are wanted (Legend see fig. 22 a) 
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Fig.. 43 a 

LEGEND for fig. .‘43 

+++ Ideatif icatioD 90% .■*. 

++ Identification ^ 90% 
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in special cases 

o sequential imagery 

S Images taken in Summer 

A Images taken in Autumn 
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CN Natural Colours 
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False Colours 
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LEVEL’ I 

LEVEL II 

LEVELIIII 

LEVEL IV 

2 . agricultural 

acres 

— grain — 

— wheat 

land 

. 

— truck crop 

f 

rice . 


plantations 

— vineyards 
orchards 




oliveyards 

citrous 

nut-plantations 

1 1 ■ 

S 

4. woods 

deciduous - 

deciduous — 

evergreen — 

— deciduous 
— ' mediterranean 



deciduous 

oak 

Macchia 


coniferous - 

spruce 

pines, stone 
pines 

larches 


5* water 

lakes 

natural • — i 

artificial 

— shallow 
deep 

+ suspension 


rivers — 
channels 


6, barren land barren rock volcanic 

sedimentary limestone 

• sandstone 
marles , 
shales 

Many categories, however, are not sharply bordered because 
of a gradual change of grey tone? i*e. a zone of mixed vege- 
tation. Even if the categories can be identified, such cases 
produce a great xmcertainly and mistake in bordering different 
land-use categories, especially if the scale is large. ; ! 
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Besides methodological problems and their solution it is ■ . 

important to calculate the cartographic accuracy. ’ 

If the geometric accuracy in the scale of 1;1’ooo’ooo tjiere 
is regarded mostly no important differences between LANDSAT- i 

1 images and topographic maps. During mapping the enlarge- , i. 

* f , 

ments squares of 1o cm x 1o cm w^re rectified by using exact !. 

‘ I * 

topographic features- identical in the images and in the 
topographic maps. 

Finally the geometric accuracy was a good as it is required 
for maps of a. scale of 1i 25o’ooo. Areas of 1 mm x 1 mm, ie.|. 
25o m X 25o m, were mapped, so that the geometric 
of topographic maps, too was obtained. 

; ‘ j|'l 

The quality of new methods are generally compared with the r. 
results of former, ’'conventional" Techniques. ■ 

Especially with respec t to woods the borders taken from 
satellite images vrere more accurate than those of existing 
land-use maps. This can be explained by this reasons: 

Those areas which were mapped frdm aerial photographs lost 
their initial accuracy by the generalization necessary for 
a scale of 1:25o‘ooo. If the map is based on field work 
the borders can only be drwan roughly. 

4 . 4i_Comgar ison__with__Automatic_Digltal^Clas sif ication 

Because of the difficulties described in the above chapters 
and because of the digital original MSS-data there must be 
asked whether an automatic digital classification would not 
be preferred to the conventional methods of photointerpre- 
tation. Only a time-cost-calculation can answer this question. 

At first, there are two important facts: 


of thb 

ORIGMAt is poor 
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'! 2 

- . A Eiiropean motmtalnous region of ca 2oo*ooo km was 

mapped . This area would be covered by 6 LAKDSAT- frames. 

- ' Level I categories v/ere mapped . They include a great 

■ lot of sub categories which had to be combined. The 
subcategories result from differences of vegetation 
and soils. 

Next we must have a look at the time needed for the con- 
ventional interpretation. 

2 months mapping 5 photogeologists 

1 month photographic ' 

j laboratory work 1 photographer 

1 month grotind truth 2 geologists 

Because of the overlapping of processes the whole map was 

done in 3»5 months. 

Furthermore, at the recent state of automatic digital clas- 
sification there must be done some premises: 

- All necessary computer programs must be established - i.e. 
programs for the statistics and programs for the correc- 
tion of grey tones, which is very important especially 

! in mountainous areas with an intensiv relief and with a 
"small scale pattern of different ■ soils. 

! I 

- The computer must be capable to combine the level I cate- 
gories from the very great nmber of subcategories often ' 

j very similar to one another. 

•Although automatic digital classification is already done 
for small European areas. and just because of the intensive 
experiments on this field, it is obvious without an exact 
calculation: 

■ { 

1; The oomputer will probably classify much quickea- than' 
than a interpreter. 
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2. The classiXication pixel by pixel will be much more exact 
. than a conventional interpretation. 

3. ;iThe costs for computer time will be much higher than the 

costs for interpreters, materials and field work. 

! 

! Furthermore the high degree of resolution is not necessary 
because areas of 8o m x 8o m can not be figured out in a 
sc8ile of 1:25o’ooo or less. ’ 

In spite of the possibilities offered by statistics and by 
computers the automatic digital classification must solve 
the same problems as the conventional photointerpretation. 

A computer needs background information as well as the human 
interpreter .But it is very difficult to code all important 
background informations in a way that a computer can work 
Jith them. 

This short discussion shows that, at least for large moun- 
tainous areas in Europe, the conventional photo interpretation 
. is and, for the near future, will be the only operational 
method,. for land-use mapping. 
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MBTHODOLOGICAL INVESTIGATIONS 
( DIGITAL ANALOG INVESTIGATIONS ) 
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Methodological investigations concerning the various 
techniques of processing multispectral data have been 
emphasized during the last two years. 

Within the scope of this activities the following 
methods have been applied: 

1) Digital treatment of MSS data including 
.*» the development of software programs 

2) Combined digital analog treatment of 
multispectral data 

3) Analog electronic data processing I ■ 

I 

^ f 

4) Fourier analytical evaluation of LANDSAT 
data by an optical computer 

A Digital Image Data Processing 

Due to the relatively late delievery of LANDSAT CCT^s we 
went into the field of digital data processing only to- 
wards the end of 1973. 

The digital treatment of multispectral data involved not 
only data processing in terms of automatic classification 
but also all data handling end preprocessing techniques 
for a more quantitative evaluation of spectral information. 

The required soft-v/are development for above mentioned 
approach has been carried out at our institution without 
adapting existing programs. 

Today the ZGF is in a position not only to have a reliable 
program package for data pre-processing but also classifi- 
cation algorithms are Implemented. 

The flow-diagram for data processing and the applied systems 
‘ , is shown in Fig. 44 > v* - 
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1 . l 2 age_data__grocessing_gackage -IMAGIN; 

1,1 Preparation of data for digital evaluation 

At the beginning of the evaluation activities no computer 
programs and arithmetic models were available to the ZGF. 
For these reasons, it was necessary to create data files 
on types which could easily be handled by computers of 
different producers. The tapes could have been read in 
the original format, but this would have required the 
creation of ASSEI^'IELER -Routines for each type of computer 
used. On the other hand evaluation programs in ASSEI'ffiLER 
seemed useless for this user-oriented application. There- 
fore a conversion routine for the data conversion from 
the original 8-blt-format to a viniversal FORTRAitexecutable 
13-Format was implemented. 

During the whole activities of digital data evaluation the 
following computer types have been used: 

- Siemens 4oo4/151 
CDC 64oo 

- CDC 17oo (for type conversion from 9 to 7 track 

and vice versa) 

- IBM 36o/9o 

- TR 44o (which will be used in the future) 

Although the utilisation of different computers for the 
evaluation activity was anything but economical, it had 
the advantage that all programs have been tested and were 
executed on either one of these. Therefore program caused 
discrepancies between different systems and compilers are 
very unlikely. In fact, discrepancies never did occur once 
the programming work terminated. 

The advantages of the FORTRAN programming, although in 
execution a little slower than ASSEMBLER routines, became 
evident once the programs are given to users, or are 
integrated to different data service system. 
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Assembler 

Assembler 

FORTRAN 

FORTRAN 

FORTRAN 

FORTRAN 

FORTRAN 

FORTRAN 

FORTRAN 



Fig. 44 Flow diagram of data processing 
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1 •'2 Check of Data Quality 

LAKDSAT CCT'S have been the only multispectral information 
in digital form. Therefore until now, all activities to- 
wards the development of a data processing software have 
necessarily been based on this data. 

At the beginning we had to face problems concerning the 
data quality. 

One principal failure which eventually occur ed were parity 
errors on the original data tapes. The maximum number of 
events that occured were about 12oo parity disturbed re- 
cords in the 234o records of tape data file. In order to 
correct these data, a parity correction routine was 
created. All parity errors up to now proved correctable. 

An analog display of either one of the multispectral bands 
showed a systematically increased intensity value for each 
sixth scan line. This effect became particular evident 
for low level intensities as for example for water bodies. 

For analog processing such as density slicing, the sixth 
scan line effect was very distvu?bing. 

For the automatic classification it did not result in the 
expected reduction of quality because of its high corre- 
lation to all spectral bands. 

Therefore no special processes for eliminating the sixth 
scan line effect were applied. 


1.3 Description of IKACIN Routines (see Apendix 1) 

The software routines developed can be separated into 4 
groups representing 

- descriptive statistics 

- data manipulation 

- classification periodicity analysis 

- support routines. 
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The full program check is listed in Apendix 1 . 
The programs involved are; 


1, Category: 

Descriptive statistics 

1.1 HDAR 

Calculation of mean, standard deviation, and 
variance for a set of multispectral data 
including graphic display 

1.2 KOFIR 

Calculation of correlation, covariance matrix 

1 . 3 REGR 

Determination of the regression fianction for 
one or several sets of data in linear, quadra- 
tic or cubic dependence and graphical display 
of the clusters and the regression function 

2. Category: 

Data Manipulation 

2.1 WIFU 

Weighting of spectral intensities for contrast 
enhancement and Illumination (for graphical 
presentation) 

2.2 KAXAN 

Principal Component Analysis for the selection 
of those axes with the high variance and the 
controlled illumination of those containing 
with the highest degree of systematic errors 

2.3 FAKTR 

Factor analysis -Principal components 

FAZEN 

Factor analysis -Centroid method 

VARMX 

Varimax rotation 

\ 

3. Category: 

Classification 

3.1 DISK 

Discriminant analysis for the supervised 
classification of multispectral signatures 

3.2 AUKLA 

Clustering for the unsupervised classification 
of multispectral signatures 


REPRODUCIBILITY OP TIJR 
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4. Category: Periodicity analysis 

4.1 FORAN Fourier analysis for the determination of one 

or two dimensional periodicities within an image 
(of only one spectral region) 


5. Category: Support routines • 

Supporting subroutines for matrix inversion, 
the generation of gray scale naps, matrix 
multiplication etc. 


On the basis of above described software package, LAKDSAT 
data were processed under the following aspects. 

1 ) Statistical evaluation with respect to the 
information content of KSS bands 

2) Data Reduction 

3) Automatic classification 
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2. Selection of Data Material and Test Sites 


Concerning OTor activities on digital data processing 
during the last year, we were mainly involved in learning 
about the significance of the spectral information gathered 
by the MSS scanner and in developing and updating software 
programs. For this basic research and software development 
the selection of test sites has not been carried out 
with special regard to our LANDSAT test area in Italy, but 
more under the following requirements . 

they should all be very close to the location of ZGF 
(not more than 15o km). Thus it would be possible to 
visit them regularly. Then only a short time would be 
needed for comparison of results obtained by the 
computer evaluation to the real situation; 

. - they should show significant featxires for the central 
European user; 

- to obtain results concerning the resolution to be 
obtained by a pixel-by-pixel selection, areas were 
chosen which contain features of high neighouring 
contrast. 

All test areas were chosen of the same size and the same 
geometry: 98 lines and 12o pixels per line = 1176o pixels 
per test area. This geometry of the frame has been selected 
simply for particular reasons of computer display, not 
for any other purpose. 

The areas in particular are as follows (fig. 45,46 ) 


Scene 

Nr. 

Name of area 

Data taken 

1. 

Ostersee Area 



2. 

Amaersee-Korth 



3. 

Hofolding Forest 

ERTS 

1o21-o938o 

4. 

Landsberg Fields 

Date: 

13 Aug. 72 

5. 

Valli dl Comacchio 

ERTS 

1218-09341 

6. 

Marina dl Carrara 

Date: 

* 26 Feb. 73 
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Fig. "IS LANDSAT-scene ID 1021-0938o with test areas 
fot the statistical evaluation 
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Fig, 4.6 LANDSAT-scene ID 1218-09341 with test areas 
for the statistical evaluation 




original 
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Scene 

• 


Area Features 


1 


Lakes, Karsh land, Coniferous forest, 
Pastures, Agricultural Areas 
Lakes, Reed Grass, Coniferous Forests, 
Leaf Tree Areas 

Agricultural -treas, Meadows, Mixed Forest 
Types 

Mixed Forest Type, Agricultural Areas 
Rice fields 
Coast area 


2 


4 

5 

6 . 


In order to look at the statistics of the nultispectral 
scarjoer data additional larger areas out of above LAMDSAT 
scenes v;ere selected (Fig. 46- ALP) 


3 * 

The statistical evaluation began "with an intensity 
(dynamic range) comparison of different data channels to 
obtain some idea of the intensity population of a tj^pical 
Central European scene, represented by the Starnberg Lake 
area. By plotting these data for two MSS charuiels in a 
Cartesian coordinate system, a distinct linear relationship 
bet 'een the intensities of ERTS infrared channels MSS 6 
and MSS 7 became evident (Fig. ). 

Assioming a straight line as the best correlation curve 
for the data set, the constants A^ and A -2 in the linear 
equation for channels 6 and 7 


were found to be 5.3o59 and I. 5 I 0 I, respectively, for this 
particular region. 


It = A I + (Aj • !■? ) 


HEPRODUCIBILITY OF 
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Fig. 47 Intensity Plot, MSS 6 versus MSS 7 


These constants were obtained bj' a least-squares method 
on the basis of the following relations; 

(2U(2;i-,’)-(2U)(^l6l7) 

" NSi.' -(ii.)' 

NIIJ, - (iU) ( 2 I 7 ) 

~ N2U-(S1,)' 


'The deviation in the measured intensity values and those 
calculated with the above formula as a function of the 
measured values of channel ? is shown in Figure 48 , The 
mean deviation v/as found to be 2,2 for a scale of 128 
intensity values for LAI'IDSAT MSS data. 
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Fig 48 Popiolation as a function of difference D 


The graphical presentation shows very good agreement with 
a nornal distribution, with populations of: 

68 in the region of 0 + 1 S 

93 9^ in the region of 0 + 2 S 

- 98 % in the region of 0 + 3 S 

with 


s =V 


“ meaMired ~ ^*ical c 



A histogran of the measured data of channel 6 and those 
derived nathenatically from the experimental values of 
channel 7 shows almost an equivalent population distri- 
bution to that presented in Figure49 for the Starnberg 
Lake are and the Italian Mediterranean scene. 


reproducibility of IH 
ORIGINAL PAGE IS POOR 


2ENTRALSTELLE FUR GEO - PHOTOGRAMMETRIE UNO FERNERKUNDUNG 


blatt 92 









n r 



TW 


h>4 



I I I ^ L_l 1_J L_1 1 
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5 Calc “ 5.3059 ♦ ( 1.5101 « l7.measuredl 


Fig. 49 Histogram of the popxilation of experimental 
and calculated values of Ig= f(ly) 

Based on these results, the analysis vras .suhsenuently 
extended to complete Central European and Italian scene.s 
indicated in Tables 1 and 2 . Evaluation of these data 

shows almost equivalent constants and A 2 for the sane 
tYpe of landscape and landforn. 


Table I. 


Scene 

Nuinber of 
pixels 

Channel 6 

1. s» 

Channel 7 
I7 S 7 

Linear equation 
A 1 Aj S 

ALP ll 

974400 

35.75 

10.83 

18.79 

6.8 

0.28 

1.84 

10.35 

ALP 31 

324800 

36.98 

13.64 

20.44 

8.72 

0.3 

1,75 


ALP 41 

97080 

44.23 

13.70 

25.19 

9.07 

3.40 

1,17 

3.71 

ALP 22 

950040 

38.21 

16.55 

20.03 

9.22 

0.3 

1.85 

26.55 

ALP 32 

950040 

40.73 

21.08 

20.62 

9.89 

0.12 

1.95 

26.55 

ALP 42 

946533 

37.71 

16.57 

20.26 

9.54 

3.54 

1.11 

10.18 

ITAL21 

936767 

22.96 

6.55 

12.18 

3.95 

0.13 

1.84 

12.1 

ITAL31 

943467 

24.12 

3.86 

11.82 

2.45 

0.11 

2.0 

17.57 

ITAL32 

943879 

25.41 

4.16 

13.08 

2.88 

0.14 

1.91 

16.19 

MED 12 

40000 





5.01 

0.0678 

5.5 


BtTBODCCBlLipf 
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Table 2. 




Scene 

Identification 


ALP ll 

Northern half of strips 1 , 3 and 4 

ALP3I 

ALP 4 1 

of ERTS 1021-09380. 13 August 1973 

ALP 22 

Southern half of strips 2, 3 and 4, 

AL? 32 
ALP 42 

ERTS 1021-09380, 13 August 1973 

ITAL21 

Northern half of strip 2 

ERTS 

ITAL31 

Nortliern half of strip 3 

1218-09341 

1TAL32 

Southern half of strip 3 

26 February 
1973 

MED 21 

Mediterranean water surface of southern 
strip of HRTS 1218-09341, 

26 February 1973 


Linearity alone, hov;ever, is not a sufficient identity 
criterion. Much more important Is the average deviation 
hetv’een the measured intensities of channel 6 and those 
derived mathematically from the experimental val.ues of 
channel 7 . The data obtained show that the average de- 
viation for the eo_uivalent landscape and land-use scenes 
are almost identical. This is valid not orUy for the Cen- 
tral European summer scene, but also for the Italian 

winter scene. The additive and multiplicative constants 
for the linear equation show a strong deviation for equi- 
valent scenes only in those cases vrhere the test area 
considered largely contains a single viniforra feature. This 
was the case for scene ALP 41 (water surface of the Chiea- 
see Lake), scene ALP 42 (mostly cloud covered), and KE3D 12 
(uniform water surface of the Mediterranean). Applying 
these results to a scene of the same type, but containing 
several different featvires, results in a ver^r wide deviation 
. betv;een the measured and mathem.aticallv derived values of 
charnel 6. 

From these results, it can apparently be concluded tha"^ a 
strong correlation, and therefore rediindancy , occurs bet- 
ween the intensity values of channels 5 and 7 for each 
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individual data point. Such, strong dependance has further 
consequences v;hich find their impact in the Tr.athercatical 
procedxires for feature recognition and identification. 

The idea vas to present the data from the four TTSS channels 
in four-dimensional mathematical space and, by identifica- 
tion of particular characteristics to recognise the presence 
of certain features in the scene. 

The fact that there is a strong correlation tetv:een the tvo 
infrared channels means that one set of data doss not 
reall^^ give nev information, so that only three channels 
and three independent data sets of the four ass'Jimed arc 
actuall3’' available. This reduction in information became 
very obvious for the factor analysis and is evident in the 
correlation matrix obtained (Table 3 ). 

Although these results could be Interpreted such tha"^ the 
chajonels 6 end 7 could be considered redundant, great 
restrictions have to be applied to this conclusions. The 
scenes considered v.'ere either vinter scenes (ita''.'') or r.i 
summer scenes (Central Europe). Both landscapes are charac- 
terised b^' The presence of highly uniform land-use natterns , 
containing features of the harvest or post-harvest period. 

To arrive at a definite conclusion, extension of this 
analysis to a greater set of land-use t^'pes over all 
seasonal periods, and particularly over the maturation 
periods of certain agricultural fruits, is necessarT*. 

It seems that in these cases ^ , as veil as for such 
featxires as vaters of different murbldity, cloud and. 
snow cover, there is no redundancy betveen the tvo infra- 
red channels, and therefore one of them cerjiot be dis- 
pensed with. 
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Table 3 Correlation coefficients LAIJDSAT MSS bands 



KSS-4 

KS3-5 

M3S-5 

MSS-7 

MSS-4 

1 . 000 

0,762 

0,213 

0,141 

MSS-5 

0,762 

1 ,ooo 

0,471 

0,333 

KSS-6 

0,218 

0,471 

1 , 000 

0,987 

MSS-7 

0,141 

0,388 

0,967 

1 , 000 

;^-ile 

working on the 

statistical 

evaluation 

of the Centra! 

Europe 

;an and Italian 

scenes, an 

interesting 

feature con- 


cerning only the water bodies becane evident, wliile the 
intensity le\’’el of European lakes ranged from 0 to 9 in 
charjiel 7, the dynanic range for the Mediterranean water 
body v.^as 0 to 2 (Fig. 5o .) 
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Fig 5o Compasrison of d3mamlc ranges 


I7 Boccq di Magro (ERIS-A 26fetrlQ73> 


\^ Sternberg Lake (ERTS-A 13 Aug '972) 


Ifi Bccca di Magro (ERTS-A 26 Febrl973) 


^6 Starnberg Lake I ERTS-A 13 Aug 1972) 
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Above described phenoniencr of the correlative property 
of MSS bar.ds led to a more detailed evaluation on xhe basis 
of special selected surface features. Furthernore this 
approach led us to discussions about the possibilities of 
laking advantage of the correlative prorerty by reducing 
the red\indant information for data compression piu'poses. 

The selection of test areas has been di'^ected at large homo- 
geneous ob.jects such as coniferous forest, v/ater 
large homogeneous vegetated areas. Furthermore, surface areas 
composed of different sub-categories of earth scientific 
interest such as agricultural field pattern, vrere also 
selected (Fig. 51 ). The graphical representation, Fig* 
illustrates the correlation of the LAITTSAT bands v;ith 
the above-mentioned surface features. 

Eased on this comparison, the follovring conclusions car. 
be reached: 

- A very high rediindancy exists betv;een the infra- 
red bands, expressed by a correlation factor 
above o ,9. 

The degree of correlation is approaching a uni- 
forn value for all samples except for the mixed 
forest. 

- The visible channels are also characterized by 

a relatively high redundancy. The absolute degree 
of correlation is more dependent on the various 
surface featxu'es. 


- Lov; redundancy exists betv.'een the visible and 
infrared channels. The variations in correlation 
are strongly Influenced by the selected sample.? - 
So minimum correlation can partly be observed 
betv/een MSS 5 and K3S 7. 
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Tl-iese examples sho’.'.- very clearly that, at least on the 
basis of LAIJDSAT, the significance of spectral bands is 
strongly dependent on the surface feature, or r.ore g'^’nc- 
rally, on the earth scientific problem in question. 

With respect to the relationship described above in general, 
it is not possible to satisfy/ different earth scientific 
requirerents v.^ith only a limited, number of spectral bands. 

The result would be the implementation of a multi s-'ect’^al 
scanner operating a high number of spectral bands. This 
requirement leads to non-man&geable data rates, esreciallv 
with respect to operational systems. On the other hand, 
the more or less high degree of redundancv between multi- 
spectral data enables the application of redundar.c 3 '’ supression 
techniques. The strong realtionship betvreen correlation 
factor and surface feature calls for a supression technique 
that is ccrried out in relation to the corresponding user 
oriented problem and that can be performed after data 
acquisition and before data transmission. 


4.1 Data Compression by Principal Components Transformation 

A well kno^^7n data compression technique is the so called 
Principal Components Transformation. This statistical method 
defines a new coordinate s?''stem describ.tng maximum dis- 
persion of the input data. The nev; coordinates or principle 
components represent a linear combination of the original 
coordinates or spectral charuiels. 

Application of this "technique on multispectral information 
indicates, that the same information, describee by the 
entire range of spectral bands, can be e^qoressed by a 
very limited number of principal components. This effect 
is •orimarilj’’ based on the correlative property of multi- 
spectral data. With respect to the number (n) of available 
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variables represented by the snectral the Irincipal 

Conponents Transformation operates in an n-dimer.sional 
space. 

From a statistical point of viev, the definition of the 
nev/ axis is based on the dispersion created by tne chance 
occurrence of N categories or data populations. Iherefors 
the dispersion or information content of a single pop'ola- 
tion and of under-represented populations of earth 
scientific interest cannot be described in an optinun ••o'' . 
Due tc this effect, data reduction on the basis of the 
Principal Components Transforcation is of necessity a 
compromise betv/een the information content of the total 
input data and the information content of minor ro’oula- 
tions. 


4.2 Application Oriented Data Compression 

The effect described above becomes crucial if the Prin;' : al 
Components Transformation is considered as an in-flight 
data compression technique for operational satellite 
systems. Furthermore as a statistical method, it requires 
elaborate computer hardv:are. 

In realit 3 ’’ an operational S 3 ’stem evidentl'^ cannot sp~lsf^' 
the entire multi-disciplinary user community', at least 
not at the same time. It will be able, ho’./ever, to cover 
a limited nur.ber of relevant applications . ^or these 
rea.sons in defining earth scientific problems, it is 
necessar 3 ’’ to have optimum information in terms of appli- 
cation oriented lovr data rates. This reouirem^ent can be 
partl3’’ fulfilled on the basis of a so-colled Supervised 
Principal Components Transformation. 

The basic idea of this method is., to_def ine_a_priori (the 
principal components of onl^' a single- or a group of 
featiires v.'hich is of interest) and to_transform_ (reduce) 
the incoming multispectral information with respect to 
these defined components. 
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3y means of this supervised reduction above-. r.entione'd 
"comprcmise-eff ect" of the Principal Conponents Trans- 
foraation is mininiaed vith respect to the disnersion of 
single or minor data populations. Furthermore, the super- 
vised transformation does not require elaborate on-board 
computer facilities. Therefore, from a technical point of 
vie'.', if offers a high potential as an in-flight data 
cacmpression technique. 


In order to demonstrate this hind of supervised transfor- 
mation, as a first step, the loadings of the orlncipai 
axis of an irohomogene agricultural area and of a homogeneous 
coniferous forest, v.^ere calculated. 

Principal Components Factor Loadings 


Agricul 

tural area: 


Conifer 

ous forest: 


I 

II 

III 

I 

IT 

MSS 4 

0.315 

-0.553 

-0.1dS 

0.898 

-0. 5 '9 

nss 5 

0.546 

-0.749 

- 0 . 1 5o 

0.927 

-0.4-7 ^ 

i-:ss 6 

O.S45 

0.528 

o.o1o 

O. A4f: 

0.511 

MSS 7 

0.754 

0.652 

0 . o24 

0.275 

o . 476 

As calculated 

above, the 

three nri 

ncipal components for 


the agricultural area contain 99 >7 'rl of the total informa- 
tion represented by the U !-;SS-bar.ds . Because of ■^he homo- 
geneous signature of the con.iferus fores'^-, ' of the 

total information can be described by only tv/o factors 
(Fig. 52 )• 
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Fig. 52 Distribution of infornation on principal com- 
ponents for agricultxare and coniferous forest 


In a second step the KSS.data of the test area (Fig. 53 ), 
characterized basically by a coniferous- and 'mixed forest 
cover and also by an agricultural area, v;as translated. 




Fig 53 Representation of the test-site (LAKDSAT-MSS 7) 
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This transformation shows very detailed inf<^rmation on 
field patterns; whereas, it is not possible to define 
the homogenity of the coniferous forest. 


A transformation of the same original data with the 
loadings of the first factor representing coniferous samp- 
les (fig. 55) was also carried out. 

This time the forest cover appears very well defined and 
the lack of information on agricultxxral areas is evident. 
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Fig. 65; Principal component image optimized for forested 


areas 


- RmbDUGBiLrrY of ths 

nmaWAL PAGE IS POOR 



ZENTRALSTELLE FOR GEO - PHOTOGRAMMETRIE UNO FERNERKUNDUNG 


BLATT l04 


By means of corresponding transformations of the original 
data through available supervised factors, the nev 7 set of 
data can be regarded as compressed input data with optimum 
information content for specific applications. 

Furthermore, by applying the remaining set of data (com- 
ponents) , a reduction of computer time in the ratio of 
9:16 could be maintained with no reduction in classifi- 
cation accuracy as the processing time req.uired is a 
quadratic function of the number of spectral bands. 


Another type of data reduction can be obtained through 
application of factor analysis (FAKTR) , in order to select 
those spectral channels with the highest degree of 
significance £>r the identif ication of certain surface 
features. Using this method for LANDSAT-data , it was found, 
that only three out of the 4 bands showed sufficient and 
linear independent information. The high degree of linear 
dependence, as pointed out already: above , between MSS 6 
and MSS 7 should be restated. 
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5. Automati^c^ ORIGINAL PAGE IS POOR 

Multispectral data require by its nature of representing 
multivariate information advanced processing techniques. 

Furthermore this data can be regarded as an optim\im pre- 
mise to apply fully automated evaluation and classification 
algorithms. . 

Qn'.the basis of data representing Italian and South German 
scenes, both supervised and xmsupervised algorithms were 
programmed and applied, 

5,1 Supervised classification 


In order to test and evaluate the supervised classification 
programs (Apendix 1) a test site near Munich (Fig. 45 
Scene 1 ) was selected. The classification was started v;ith 
ob;jects easily identifiable by their areal extension and ■ 
spectral characteristics such as water bodies, coniferous 
forest, \irban areas, and pastures. The result of such a 
supervised classification is demonstrated by fig.^ 56 

In order to prove the reliability of above classification 

i 

intensive groxmd checks and comparative studies based on 
aerial photographs were carried out. 

From this it turned out, that some classes which are uniform 

from, an earth scientific point of viev/ could not be fully 
recognized by the automatic classifier. Marsh land for 
example, v^hich is composed by different vegetated and non 
vegetated surface ' featiu'es could not be separated. Training 
the computer with data, representing the unhomogeneous category 
"marsh” interference with other surface features such as 
coniferous forest occured. 

Dividing the cluster of ' msishland into more subclusters it 
has not been possible to map the true areal extension of 
the marsh itself in spite of a high statistical accm'acy in 
terms of separability. Such difficulties in defining special 
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surface features are not caused by the algorithms used. 
Rather the reason for classification errors is the signi- 
ficance of the available spectral information for solving 
a partictilar lash. 

Further difficulties arise especially with respect to the 
resolution capabilities of the MS-scanner and the typical 
pattern of agricultural used areas in Europe. 


In contrast to the extensive land use, characteristic of 
North America, we are faced with an intensive management of 
our agricultural areas. This has resulted in very small 
and ’Wiomogeneous field pattern. 

Due to this situation, it is necessary to make a very careful 
selection of the picture elements for training the computer 
with certain test features. A further problem is, that the 
test site very often is much smaller than the dimension of 
a resolution element leading to a mixture of different 
spectral significances for one cluster. Furthermore this 
results in a very poor separability crucial for automatic 
classification pvirposes. 

The same problem of gathering unhotnogeneous spectral information 
occurs at the borders ' of tv^o significant surface features. 

The degree of misclassif ication for this case can be mini- 
mized by comparing neighbouring picture elements. 


Due to above mentioned problem in classifying phenomena 
"below the resolution" other approaches of automatic in- 
formation extraction techniq.ues should be taken into con- 
sideration. 

To carry out a spectral- statistical inventory of agricultural 
used areas over one vegetation cycle in combination with 
detailed ground truth reports may result in correlation which 
allow a further use of Landsat data in terms of a pure 
momerical evaluation. A premise for this approach has to 
be a maximum repetitive LArJDSAT coverage. 
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5.2 Non supervised Classification 

For carrying out a non supervised classification tvm routines 
were implemented. 

The first routine stores all data necessary for an automatic 
classification in core. As the nvimbers of storage elements 
necessary for the classification depend on the number of 
picture elements by n (n - l)/2, its execution requires; much . 
storage, but, on the other hand, is very fast. 

The second routine has in core only the information of the 
clusters found, and the information of the actual data -points. 
It is as fast as the first routine, but does not need as 
much storage, 

V/ith these routines a separation of "those clusters vith 
a very high degree of agreement has been obtained. It soon 
became evident, however, that it was not meaningful to use 
all data elements for the classification, because of the. 
great computer time necessary for these procediores. Therefore 
the decision was made to use only a very small number of 
preselected pixels of a specific significance- or to select 
the data points from a certain land-scape t 3 T?e by i-andom 
number access. 

For classification purposes the Laguna die Comachio 
south of the Po-delta was chesen (Pig. .^6, Scene 1) 

The area covers water bodies, large agricultural fields 
with xaniform plant soil association characterized by ex- 
treme changes in soil moisture and shallow water in 
irrigated rice cultures. 

By applying AUKLA and AUKLB it v;as possible to separate 
four typical c-.-Cegories. The computer print out of the 
evaluated area (Fig. 57 ) shows the following categories: 
deep v/ater, shallow water, channels, water covered fields 
or areas with extremely high soil moisture, fields with 
fresh vegetation, dry area v/ith vegetation. 
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Fig, 57 
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reproducibility of the 

5.3 The "Slope Method” ORIGINAL PAGE IS POOR 

V/ith respect to the computer time consuming classification 
operations the possible information, gain has to be care- 
fully examined, » Even .-or a purs scientific data handl.tng 
computer time often becon.es a critical factor which calls 
for operating ma^'- be less complex processing softvmre or 
analog eouipraent. 

A very easy- and still effectful possibilitj’’ of digital 
evaluation of multi spectral data is, taking into account 
the reflection properties of some dominant surface featu- 
res v;ithin the spectral range covered by the ERTS-MSS 
scanner. 


Sons basic spectral properties are the reflection capabi- 
7-ities of plants and water in the near region of the infra- 
red spectrxjm. In general, reflection characteristics depend 
on the spectral range or wavelength of the surface considered. 
For the particular case of the LANDSAT-datd -the slope has been 
defined as ■ 

a = f (contribution of visible ~ contribution of 
infrared part of the spectrum) 

L + Is' _ U G 


This technique is not based on a statistical approach. Each 
pictiore point, represented by 4 "Variables - the 4 spectral 
measurements - is classified with respect to say an certain 
angle intervalls, and assigned with a corresponding symbol 
and printed into the maps. This simple decision making 
program, does not require large amounts of computer time, and storag 
and can be regarded as a kind of digital quick look operation 
(fig. 58 ). 

The idea of the slope technique is not to imitate. density 
slicing techniques but to have the possibility to perform 
a certain kind of preselection of s^n’face features. 
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Fig. 58 Spectral slope method 

Under this aspect, it is possible to get a 'distribution 
display of different v/ater bodies but no idea about i is 
relative degree of pollution. 



Fig, 59 Color composite re- 
, presenting 3 classes 

of turbid water 
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In Fig. 59 the Versilian coast is evaluated, using this 
method. Three different classes of turbid v;ater bodies are 
distinguished. Main pollution surces SE of the harbour 
basin of Marina dl Carrara, the distribution, transporta- 
tion and degree of pollution in Thyrrenian Sea can be re- 
corded in relation to different slope angles. Decreasing 
slope is a function of the reflection and absorption pro- 
perties within the visible and near infrared portion of 
the spectrum. 
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B Analog/Pig'i'cal Treatment of ilultl spectral Data 

The results of fully automated classification procedures . 
are data in maplike form which consists of averaged infor- 
mation with respect to the' information content of the ori- 
ginal data. By this the user is confronted with a final 
data product \^hich does not permit anymore any further inter- 
pretations. 

Taking into account ^ the interference and complexity ■ of 
ear cli-scientif ic phenomena and the difficulties described 
above for many applications the researchers know how is 
indispensible . • 

Especially in an experimental phase it is of utmost inport- 
ance of being able to evaluate multispectral data by combining 
the advantages of statistical techniques and- the interpreters 
know how. 

1 

From an earth scientific point of view a possible solution 
to take advantage of both can be achieved via the follovring 
steps. 

a) statistical enhancement of class separability 

b) combination ox transformed data by operating in color space 

c) interpretation 

By this technique dev'-xsion making algorithms which actually 
result in an averaging of da" a .Is left to the interpreters 
interaction. 


a) Enhancement of class seperabilitv 


This preprocessing technique represents a necessary first 
evaluation step for both the fully automated classification 
and the combined digital and analog treatment of multispectral 


data. 
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The separatilitv of surface classes represented t;- 
clusters in an N dimensional space is comrronlv increas 'd 
by applying the principal component transformation. Since 
a linear transformation of all data points considers the 
weighted intensities of the oz'*ifinal bands, nev.’ dimensions 
or components can be determined. Due to the enhanced 
separability of the clusters, the pictural representation' 
of the new components offers enhanced innximzm contrast. , 

In order to obtain principal component images which sl.ow 
maxim.’tm information for selected phenomena or applications, 
the transformation must be optimized. The optimization is 
earn" e'* out by applying this technique on selected treining 
areas ^as it is demonstrated on page ). 

The image of fig. 61 represents aveighted transformation of 
the original bands (fig. 6o ) . It is evident that none of 
the original bands (MSS 4,5,6 and 7) have a comperable con- 
^rast in the region of field patterns. 

'Ith z'espect to the ''■arious application-dependent •^ecr’ire- 
merts, it should be possible to define standard transiorma- 
tior by taking into account also changing envirenmenta?;. 
parsmeters. "herefore this approach has to be based on 
along term exper:* mental remote sensing mieasurement program. 



Fig.ao LANDSAT-MSS 4 




«TnTlJTT .TTY 






Fig. 61 Transformed image (first principal component) 
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The tiansformed spectral information now can be used as a 
basis for a further ^utona-^fc classif ication. With respect 
to a further analog treatment, a selecticn of three trans- 
foriuet images (either the first three principal components 
or a selection of optimized transf ormetions) cr.n be u.' 
for the generation of color joirposites. '/.'ith r'espect to JiC 
G>roanied information content of the black and vh.ite incut 
i;nn.5:es the resulting perceptable color ■fifiarerc&s are ve’^y 
Indicrtive of the distribution of charajteristic s'oi'fect: 
signatures. 

For ceconstretion some areas cut of an FFTS-fre" . shc-ing the 
Alps brtv’een I'xinich and Bozsi ••ere selected '"fif. 62 , liarel 

and iL ). The image of fir. 62 illustrates the nlarred regivT. 
of the rrurnauer Moos. The image is composed from KSS band t 
(bine' , 5 (green) and* 7 (red). ‘This false color imag.': shv. a 
a rele tivel’’’ low degree uf different: aticn. Tv•^ t.ai.n 
gories, grassland (light red) and bog areas can be seperated.. 
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Fig. 63 False color composite "Mumauer Moos" 


The color composite (fig. 64 ) derived from statistically 
tz’ar.sformed data is characterized by a much higher contrast 
and an enlarge^ “spectrum" of colors. 

Plant ecolopist.s could derivi from the tra isf ormed image 
very good inf o.'-’ otion on the a eel diat^lbuticn cf raised- 
tog (blueish) and trasition bog A differenclation 

v’hlch could not be obtained on ttie falss-c.'lor presentation. 
Furthermore tvc T,.'pes of i, Itgrr t 'ov*' and green can 

be seperated. 
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Fig «64 Color transformation on the basis of statistically 

preprocessed data 

Lue to the moment an eous experimental status, it ’vould he 
desirable to quantify and standardize color coding of opti- 
mized information on the basis of 3 dimensional color space. 
The term standardization means a transformation inxo per- 
ceptable color differences. 

Quantification means a renroducabie relationship between a 
distinct phenomena and color. Such relationships are for 
example red color toner, indicating veneration if false color 
images are interpreted. 

Such "quantified data" can he obtained if the input 
data are derived from statistically optimized preprocesses. 
For demonstration the data of test area 2 (Figi>2 ) 
were selected. 

The following color transformations show enhanced infor- 
mation due to different- svirface phenomena (Fig. 65 ,66 ,67 ) 
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Fig. 68 Aerial photograph exhibiting the same scene 
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On the basis of above processed information the user would be 
provided with optimum data. In comparison to the computerized 
classification techniq.ue, the transformation into the color 
space procedure, results in data still "untouched” by decision 
making errors or averaging effects. 

The decision making aspect still left to the interpreter 
coiald be carried out via isochrome mapping techniques or 
even by an interpretation based on the earth scientific 
know ho\f. 


UE3PRODUCIBILITY OF THE 


3. Analog Electronic Data Evaluation 

The increased interest of scientists in Landsat data and 
their requirements towards easy to operate and interactive 
data processing possibilities led to the implementation' of 
an analog data processing system (Phase I, completet July 
1975 ) t v/hich will be interfaced to a minicomputer (Phase 
II, completed beginning of 1976). 

By means of the interface combination, the system offers 
hybrid capabilities. For quick look data evaluation based 
on photographic LANDSAT-data~products it is possible to 
operate the analog electronic circuits seperately. The 
System has the following standard evaluation f\mctions: 

I Analog Electronics 

- point-vrise density reading 

- densitjF slicing 

- line by line intensity level analysis of the 
whole image including selection of particular 
features 

- combination of multispectral and 'sequential data 
with respect to follov/^ing algorithms 

ratios; products, sums, differences 
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- color space operations with original or preprocessed 
data 

- electronic representation of scatter diagrams of 
any tv/o data sets. 

The digital operating circuits will he supported by problem 
oriented software routines covering the following function, 

~ CCT input capabilities 
~ geometric correction of LANCSAT data 

- classification statistics 

- basic image enhancement. 

The results obtained by a pure digital treatment of data can 
be displayed on a line printer or via the analog electronic 
circuits. 



Fig. 69 Scheme of an analog/digital data processing 
system 
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C Fotirler analytical evaluation of LAMD5AT data “by 
an optical com'puter 

Diiring the last year coherent optical processing tech- 
niques were applied on LAM)SAT-data under various aspects* 
The optical bench used is available at the Institut ftir 
Hachrichtentechnik at the Technical University of Munich. 

Optical processing was carried out for the following 

f 

purposes: 

- image manipulation 

- quantitative measurement of spectra 

Image manipulations were done with an experimental system! 
exhibited by Pig,'7o 



Fig,7o Schematic representation of optical bench 
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Using a coherent light source and a corresponding lens - 
system, the two-dimensional pattern (Image) of the input 
plane is transformed into the Fourier transform (fre- 
quency plane). Now it is possible to remove spatial fre- 
quencies from the Fourier spectrum by applying optical 
filters. 

Afterwards the filtered image can be optically rebuilt 
(output plane) which again requires a corresponding lens. 

By applying high- or low pass filters image enhancement in 
terms of contrast manipulations can be carried out. 

More intersting for geological applications is the so called 
directioneal filtering. 

For this pircposes a 10° wedge filterv/as used to remove 
those directions which are easily detectable by the eye 
and which may camouflage minor directions. 

By using an S 190 A enlargement representing the Apenine 
foreland between Pascara and Ancona (Pig. 71 a) , an 
attempt was made to remove especially those directions 
that are not necessarily of tectonic origin in terms of 
fracture control. The orientation indicated is due pri-$ 
marily to the drainage system direeted perpendicular to 
the Adriatic coast (Fig.?ia). The reconstruction (Pig, 7iij) 
shows very distinctly the prolongations of the main llne- 
ations which cut the Apenines in an ENE and anaESE direc- 
tion (Fig^. 7le) .Furthermore, a secondary system acompany- 
ing the main structures can be detected, especially in the 
northern part of the area selected. 
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Eastern Apennine foreland 
betv/een Ancona and 
Pescara 

SKYLAB S J30 A, Sept. 73 
( Beta i 1 ) 

Reconstruction from an 
optical Fourier trans- 
f o r;ii 


( a ) 


. Pescara 





( b ) 
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Of special interest for tectonical applications is to ob- 
tain a quantitative measure of the directional distri- 
bution of linear features. 

The conventional way of obtaining such information re- 
quires time consuming evaluations and measurements of length 
and direction of each single lineament. The results are 
then statistically plotted in form of rose diagrams. 

In order to evaluate the two-dimensional spectrum in a 
quantitative way, the relative brightness of the spectrum 
has to be measured. 

The information of brightness - distribution was obtained 
by using a photo multiplier which was connected to the 
frequency plane. By a step wise rotation of a narrow wedge 
filter the intensities can be easily recorded with respect 
to the direction. 

In order to obtain information on the accuracy of this tech- 
nique a lineament (Fig. 72) was optically evaluated. 

The result was compared with a conventional evaluation. 
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Comparison optical/hand evaluation of a li- 
' ^ neament map 
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In a next step a small part out of a LANDSAT scene (Fig.74 ) 
was used. The distribution of image frequencies were di- 
rectly measxired and again compared v;ith the lineament- 
distribution obtained by a pure interpretation. (Fig. ) , 



Fig. 74 LANDSAT MSS-6 Rhine Graben with test site 
for optical processing 
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Fig. .75 Comparison optical/hand evaluation of testssite 


This figure demonstrates again, that it is possible to 
obtain information at least on the distribution of domi- 
nant linear features by applying optical processing tech- 
niques. 

Above described good possibilities led to a more detailed 
evaluation taking into account the frequency distributions 
which can be derived from different spectral bands. 

Fig. 76 demonstrates the frequency distributions of the 
same test area (Pig.74 )v/ith respect to MSS bands 4,5,6 a. 7. 
A good comparisons af above readings is given by the 
characteristic peaks originated from the scan-lines. An 
effect which is distinctly expressed on MSS 4. 
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Fig. 77 Comparison Intensities band 7 versus band 4 


These two characteristic intervals could be related to the 
dominant fractxrre system in the vicinity of the Rhine Graben. 
Category I corresponds to the Rhine Graben-parallelfractxire 
system and Category II to the lineations which are + per- 
pendiciolar to the parallel structires. 

This example has to be regarded as an attempt to com- 
bine nmltispectral- with spatial-? information on a more 
quantitative way. 
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